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The  research  which  was  conducted  under  this  contract  fell  mainly  into 
three  categories:  gain  switching,  injection  locking,  and  photon  statistics. 
In  addition,  some  work  in  the  foundations  of  radiative  transport  theory 
and  spectroscopy  was  also  performed.  A  brief  summary  of  this  work  will 
now  he  given. 

I.  Gain  Switching 

The  use  of  a  metastable  state  as  the  upper  laser  level  offers  the 
possibility  of  obtaining  large  inversions  and  thereby  storing  large  amounts 
of  energy  in  the  lasing  medium.  The  problem  associated  with  this  is  that 
the  linear  gain  which  one  obtains  from  the  transition  from  this  level 
to  the  lower  laser  level  is  quite  small.  The  possibility  of  "dumping" 
the  energy  stored  in  the  mctastable  level  by  the  application  of  an  external 
field  was  investigated.*  The  external  field,  by  inducing  a  dipole  moment, 
increases  the  transition  probability  between  the  two  levels.  One  can, 
therefore,  obtain  a  large  inversion  and  then  switch  on  the  field  to  obtain 
a  good  value  for  the  linear  gain. 

An  estimate  of  Ml,  h2,  and  induced  Cl  transitions  (for  an  applied 
electric  field  of  10~’  Vcm  !j  indicated  that  transitions  satisfying  the 
following  two  requirements  would  be  good  candidates  for  gain  switching: 

A)  they  are  forbidden  for  Ml  radiation 

B)  they  occur  at  wavelengths  of  1  to  10  ..m 
Vibrational  transitions  in  homonucleur  diatomic  molecules  satisfy  these 
conditions,  lor  I!.,  one  can  obtain  a  gain  coefficient,  K,  of 

..  AN  -2.  -1 

k  =  —  (5  x  10  )  cm 

where  AN/\  is  the  relative  population  inversion  and  the  electric  field  is 
a  imed  to  be  10'’  V/cm. 
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If  one  tries  to  use  electric  fields  much  larger  than  those  considered 
above,  breakdown  will  occur.  For  the  alternating  electric  field  produced 
by  a  giant  laser  pulse  it  is  possible  to  obtain  much  higher  values  of  the 
field.  This  would  allow  the  gain  to  be  increased  by  two  orders  of  magnitude 

Spin  forbidden  transitions  were  also  examined  to  see  if  they  would  be 
of  use  in  gain  switching  schemes.  In  particular  che  transition  from  the 
state  of  molecular  oxygen  to  the  ground  state  was  considered  because 
this  state  can  be  produced  in  chemical  reactions  so  that  a  practical 
pumping  scheme  is  available.  An  inhomogeneous  magnetic  field  will  cause 
spin  multiplet  intercombination  and,  therefore,  make  spin  forbidden  transitu 
possible.  Actual  calculations  for  the  ‘A  -  oxygon  system,  however,  showed 
that  the  necessary  magnetic  fields  were  not  practically  obtainable. 

II.  Injection  Locking 

There  are  applications  where  it  is  necessary  to  have  a  high  energy 
laser  with  a  particular  polari cation .  Due  to  practical  considerations 
it  is  not  always  possible  to  use  intracavity  polarication  selective  elements 
Die  laser  output  is  then  unpolar i cated ,  or  because  of  resonator  geometry' 
or  reflective  properties  of  the  intracavity  optical  surfaces,  the  laser 
may  have  a  prefered  polarication  which  is  usually  not  known  apriori ,  and 
in  some  cases,  is  time  varying. 

It  is  possible  that  injection  locking,  which  has  been  successfully 
applied  to  frequency  control  lasers,  may  also  be  used  for  polarication 
control.  In  fact,  experiments  have  been  performed  where  a  polarized 
output  was  obtained  by  injection  locking  an  unpoluriced  rare  gas  iiulide 
laser.  We  wanted  to  determine  if  polarication  selection  by  injection 
locki  will  work  for  a  wider  range  of  lasers. 


We  have  extended  the  seiniclassicul  laser  theory  to  include  laser 
operation  with  different  polarizations  and  with  an  injected  signal.  We 
considered  two  problems.  The  first  is  the  effect  of  a  polarized  injected 
signal  on  an  unpolarized  laser.  The  second  is  the  effect  of  a  polarized 
injected  signal  on  a  laser  that  has  a  prefered  polarization,  he  wanted  to 
know  if  the  injected  signal  can  force  the  laser  into  operating  with  a  different 
polarization.  This  latter  case  will  be  of  interest  in  the  scaling  of  lasers 
by  using  a  multiplex  array.  It  may  also  find  applications  in  lasers  with 
annular  gain  regions  which  at  present  operate  with  axicon  type  resonators 
that  produce  outputs  that  have  time  varying  polarizations.  Our  goal  was 
to  obtain  as  understanding  of  the  mechanisms  underlying  polarization  control 
via  injection  locking.  As  such,  we  have  not  dealt  with  tue  more  computa¬ 
tionally  complicated  though  possibly  more  useful  cases  of  operation  in  the 
strong-signal  regime  and  with  an  unstable  resonator.  We  are,  at  present, 
investigating  the  above  two  problems. 

Our  results  indicate  that,  except  when  the  injected  signal  and  free 
running  laser  polarizations  are  orthogonal  to  each  other,  it  is  not  possible 
to  rotate  the  laser  polarization  so  that  it  lines  up  with  the  injected 
signal  polarization.  However,  the  difference  in  polarization  may  be  made 
arbitrarily  small  by  increasing  the  injected  signal  strength.  The  reason 
for  0  ^  Oj  for  finite  injected  signal  strength  and  0<0.<n/2  may  be  seen  b. 
resolving  the  injected  signal  electric  field  into  a  component  that  is 
parallel  to  the  free  running  la.scr  polarization  and  component  that  is 
perpendicular  to  it.  The  parallel  component  has  a  greater  effect  per 
electric  field  strength  than  tire  perpendicular  component  because  of  the 
higher  laser  gain  in  that  direction.  Consequently,  the  final  laser 
polarization  will  always  he  between  the  free  running  laser  and  the  injected 
signal  polarizations.  When  0  ^  =  '■'/ 2 ,  there  is  no  parallel  component  ot 
the  injected  original  and  therefore  it  is  possible  to  make  u  =  6.. 


If  we  consider  the  situation  where  we  have  an  array  of  lasers  and 
where  each  laser  has  a  different  preferred  polarisation  (i.e.  the  output 
of  the  free  running  laser  array  is  unpolarized),  then  the  distribution 
of  polarization  directions  of  the  phase  locked  array  w i 1 1  have  a  width 


given  by 


A0 


_  1_  A»  /  l  ” 
2  Ac  \  4 


where  Act  is  the  difference  in  the  gain  for  the  components  of  the  polarizati a. 
parallel  and  perpendicular  to  the  injected  polarization,  Ac  =  c/2L  where 
L  is  the  length  of  the  laser,  I  is  the  intensity  of  the  laser  inside 
the  cavity,  and  4  is  the  intensity  of  the  injected  signal  outside  the  cavi,./ 

III.  Photon  Statistics 

Our  group  has  had  an  interest  in  photon  statistics  problems  for  sev.-  al 

/ears.  Recent  work  includes  a  generalization  oi'  /rev Lous  results  for  a 

single  mode  m-photon  absorbtion  process  to  an  arbitrary  number  of  modes.- 

Two-photon  lasers  were  also  studied.  The  photon  statistics  of  such  a 

laser  with  a  loss  mechanism  simulated  by  two-level  atoms  which  absorb 

a  single  photon  was  found.-5  This  was  then  used  to  investigate  nonclassicai 

4 

effects  in  the  light  emitted.  Two-mode  lasers  were  investigated  as  well 
and  the  intensity  correlations  for  a  laser  of  this  type  with  coupled 
transitions  were  found.-5 

IV.  Other  Work 

fundamental  work  in  spectroscopy  was  also  done.  It  was  shown  that 
time-delayed  measurements  allow  one  to  obtain  accuracy  greater  than  the 
natural  linewidth  in  the  measurement  of  spectral  lines.6  The  foundations 
of  the  theory  of  radiative  energy  transfer  were  also  the  obj  t  of  a  stiu...  .  ' 


Wave  propagation  in  a  random  medium  was  considered  and  various  quantities 
appearing  in  the  usual  phenomenological  theory  were  related  to  the  stochastic 
properties  of  the  medium. 
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We  consider  situations  whereby  energy  could  be  stored  in  a  metastable  state  and  then  "dumped"  by  switching  on  an  ex¬ 
ternal  1  ield  thus  enhancing  the  coupling  to  a  lower  state  at  some  later  time.  The  present  calculations  indicate  that  tins  con- 
cvpt  merits  experimental  study. 


In  many  high  power  lasers,  one  often  seeks  a 
population  inversion  between  levels  which  are  only 
weakly  coupled  in  order  that  a  substantial  inversion 
may  be  obtained  without  excessive  decay  to  the 
groundstate  via  spontaneous  emission.  In  some  cases, 
however,  this  leads  to  the  complication  that  the  linear 
gam  associated  with  the  transition  between  these  two 
levels  is  very  small.  It  is  thus  interesting  to  consider 
situations  whereby  eneigy  could  be  stored  in  a  meta- 
stable  state  and  then  "dumped"  by  switching  on  or 
enhancing  the  coupling  to  a  lower  state  at  some  later 
time.  This  might  be  accomplished,  for  example,  by 
applying  a  strong  electric  field  to  a  metastable  state. 
This  "field  induced"  effect  is  well  understood  theo¬ 
retically  [  1 ,2|.  and  has  been  verified  in  numerous 
experiments  [5 .6. 14.1  (>| .  Thus  we  propose  to  use 
field  induced  transitions  as  a  technique  for  switching 
the  gain  in  an  inverted  medium  by  the  application 
of  an  external  field.  This  technique  could  also  be  ap¬ 
plied  to  amplifier  systems.  Eurtherinore.  it  would  al- 

Resiarch  supported  m  p.ul  In  I  lie  US.  Air  Tone  Office 
o’’  Scientific  Research  under  Grant  No.  A  I  OSR-KIMI27K. 
and  In  the  M. ix-l’l. mck-Gesellseli.il t  rut  borderline  der 
U iseenscll.il ten,  Mtinchen.  West-Germ. ins  .  One  ot  us  ( AWO) 
was  supported  by  the  NST  Materials  Research  Laboratory 
Program. 

*  Senior  Humboldt  bellow. 


low  one  io  precisely  control  the  onset  of  lasing  action. 

In  this  paper  we  shall  demonstrate  the  feasibility 
of  this  concept  by  studying  a  specific  type  of  molec¬ 
ular  transition,  and  show  that  gains  of  several  per¬ 
cent  per  cm  can  be  obtained.  The  present  calculations 
indicate  that  the  concept  of  gain  switching  deserves 
further  experimental  study. 

In  the  present  calculation  we  wish  to  compare  the 
electric  dipole  transitions  which  are  induced  via  an  ex¬ 
ternal  field  (induced  dipole  IE1 )  with  other  types  of 
allowed  transitions,  leading  to  magnetic  dipole  (Ml) 
and  electric  quadrttpole  (L2)  radiation.  The  orders  of 
magnitude  involved  in  Ihese  processes  arc  well  known 
in  the  M  I  and  LI  cases  and  can  be  easily  estimated  in 
the  Ili I  case.  Lor  11:1  transitions  the  dipole  moment 
induced  by  external  field  /■'  is  given  bv  [A] 

(//.•  =  a/’’, 

a  =  S  ( b I /->  1  r > ( r | /J | a >  ( ----- )  =  \  b|a|  a> 
r  r  ~  ,:b  '  r  '  j 

(  I  I 

where  I)  is  the  electric  dipole  operator  ,  a.  h  and  t  label 
the  lasing  and  intermediate  states  in  the  transition  con¬ 
sidered.  Tensor  a  is  the  matrix  element  of  some  scalier- 
ing  tensor  operator  at  (4 1 .  Actually  formula  ( 1 )  gives 
the  symmetric  part  of  a  only,  which  can  be  shown  to 
dominate  the  antisymmetric  part  in  most  cases  The 
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scattering  tensor  may  be  estimated  in  terms  of  an 
average  transition  energy  A/:'  and  typical  dipole  matrix 
element  J  as  follows 


_  ,  <b|/)-/3|a> 
a  --  _y. 


,  ~  2  <"b‘l  I'  ~  -  ~\T-  ~  4  X  ,  o 


a/-: 


(2) 

(3) 


Here  we  have  assumed  an  external  electric  field  /•' 

=  I  O'2*  V  cm ' 1 ,  a  typical  transition  energy  A/:'  —  4 
X  10'  '*  ere  corresponding  to  radiation  at  optical 
wavelength  \  =  0.5  gm  and  a  dipole  matrix  element 
of  the  order  of  magnitude  J  ~  c</().  Estimates  ol  ihis 
type  lead  to  rigorous  bounds  on  the  scattering  tensor 
in  all  cases  where  the  transition  matrix  element 
<h|«la>  can  be  calculated  in  terms  of  diagonal  matrix 
elements  |4|  This  occurs  for  example  for  atomic 
transitions  between  fine  splitting  components  in  LS- 
coupling  and  tor  vibrational  transitions  in  molecules. 
q  max  be  appreciably  larger  than  estimated  in  (5 1  if 
there  are  strongly  resonant  levels  r  close  to  any  of  the 
levels  a  or  b.  In  genera)  however  (3)  gives  the  correct 
order  of  magnitude  for  atomic  systems. 

For  Ml  transitions  the  magnetic  dipole  moment  is 
a  multiple  of  the  Bohr  magneton  |3| : 

i/^l  |  —  cliilmc  =  \(e~  ‘hc)cau  ~  3  (i  X  10'  *  J.  (4) 

Finally  for  F2  transitions  the  rate  is  obtained  as  fol¬ 
lows  (3) :  Neglecting  numerical  factors  of  order  unity 
we  replace  the  dipole  matrix  element  J  by  the  electric 
quadrupole  moment  Q  and  divide  by  the  wavelength 
X  of  the  radiation  emitted.  Since  (J  equals  roughly  a 
typical  intra-atomic  distance  times  the  dipole  moment 
we  obtain  for  the  corresponding  matrix  element: 

J,  :  =  X''(/ =  X”l<?(,u,~  10  4  c/  for  X  =  0  5  quit 5) 

In  a  working  gain  switched  laser,  application  ol  the  ex¬ 
ternal  field  should  increase  (he  gam  coefficient  ol  the 
medium  at  least  by  an  order  of  magnitude.  It  is  clear 
from  the  estimates  given  above  that  xve  should  look 
for  transitions  which  meet  the  following  requirements, 
a)  they  are  forbidden  for  M2  radiation  (eq.  (40; 
h)  they  occur  at  wavelengths  of  I  to  10  /am,  (eqs. 
(31,(50. 

Increasing  the  wavelength  by  an  order  of  magnitude 
will  decrease  the  rate  for  F2  rclati  >  1. 1  radiation 
by  iwo  orders  and  increase  the  rate  ol  II  1  telaiive 
to  1.2  radiation  by  the  same  factor.  These  require- 


304 


merits  are  met  by  vibrational  transitions  in  lioino- 
nuelear  diatomic  molecules. 

Homonuclear  molecules  like  Hi  do  not  possess  a 
transition  dipoie  moment  for  either  vibrational  or  rota¬ 
tional  transitions.  Furthermore  their  magnetic  dipole 
moment  does  not  depend  on  the  mtcrnuclear  separa¬ 
tion;  hence  vibrational  transitions  cannot  occur  for 
Ml -radiation  either.  These  transitions  are  observed 
however  as  electric  quadrupole  transitions,  in  Raman 
scattering  and  as  field  induced  dipole  transitions  f 5 .6 J . 
Vibrational  transition  with  Au  =  1 1  are  strongest  and 
are  the  only  ones  observed  so  far.  Transition  mients 
for  |  Au  J  >  1  are  found  to  be  at  least  one  order  of  mag¬ 
nitude  smaller  [7). 

The  rotational  selection  rules  obeyed  are 

A /  =  -2  O-branch. 

A /=  0  0-branch. 

A/  =  +2  S-branch. 


For  Hi  the  transition  between  the  first  vibrational 
state  and  the  groundstate  occurs  at  a  wavelength  of 
X  =  2.4  pm.  Hence  we  can  expect  the  field  induced 
transition  rate  to  be  at  least  one  order  of  magnitude 
larger  than  the  quadrupole  rale. 

As  motivated  above  we  now  turn  to  a  calculation 
of  the  field  induced  rale  (2).  Lei  us  first  recall  the 
physical  mechanism  of  Field  induced  transition  in 
homonuclear  molecules.  The  external  field  will  in¬ 
duce  a  dipole  moment  in  the  electronic  charge  cloud. 
This  induced  moment  will  depend  on  the  internuclear 
separation,  and  is  hence  coupled  to  the  vibrational 
motion  of  the  molecule.  The  rotational  selection  rules 
stem  from  the  fact  that  in  Ihc  transition  one  quantum 
of  angular  momentum  is  transferred  to  the  photon 
emitted  and  one  quantum  is  exchanged  with  the  ex¬ 
ternal  field. 

The  square  matrix  element  of  the  induced  i;-:;.-.!- 
tion  moment  then  is  given  by  (o| : 


g77-r  £  KuAU|D|o7Vlf>|’ 

2d  +  1  MM 


+ 


4_  _  JjJ  +  lj _  : 

45  (27  +  3li  2>  -  2)  7,11 


/  - 


/  =  ./. 


-I—  7(7  t  ) )  1  1 

=  15(27  +  1X27  -  I)71'1  ‘ 


2  (7  +  2X7  +  U  1 
15(27  +  1X27  +  r>  7,>i 


7=7  +  2 
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Table  l 

l.xpenmenial  and  theorencal  results  lor  the  polarisabilus  matrix  elements  o!  the  molecules  Ms  and  N2 


possess  a 

tonal  or  rota- 
elic  dipole 
ear  separa- 
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ongest  and 
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!  will  in- 
arge  cloud, 
nternticlear 
'rational 
action  rules 
■te  quantum 
:  photon 
tli  the  ex- 

ed  transi- 


/  = 

J'  =J  -  2 

f  =j  + 


It;  Nj 


exp. 

thcv>r. 

exp. 

thcor. 

otQi  <  10~25  cm^) 

1.2  |A|. 

11  1 15  | 

0.5 

1M. 

0.48.0.38  (P 

>01  <  1 0_2S  cm3) 

0.7  1 1 4  | 

0.9  |1S  | 

0.7; 

it'!. 

0.54,  0.38  1 1' 

\_1  MO3  cm"* ) 

4.16 

2.36 

Here  aM|  and  7m  arc  matrix  elements  of  the  electronic 
groundstate  polarisabililies  a(R)  and  ->( r  1  iR  -  inter- 
nuclear  distance)  taken  between  the  vibrational  levels 
considered,  a  and  7  are  defined  by 


The  gam  coefficient  K  is  then  given  by 


k  =  AV  A: 

4-: 


W 

cM\  ‘‘I 


Q=jai+TQ»-  7  =  a:: 
where 

“  Qy  y  ~  ^  t 

arc  the  components  of  the  electronic  polarisability 
tensor.  Transitions  in  the  Q-branch  (A !  =  0)  are 
strongest.  Neglecting  contributions  from  7(]|  which 
are  of  the  same  order  of  magnitude  as  probable  errors 
in  a„|  we  obtain  for  O-transitions  the  spontaneous 
rates 


=  0.9  X  10  4  s'1 

for  1 1 ■. 

=  3.4  X  10“6  s  1 

for  N, 

both  at  /•'  =  10'S  V  cm. 

Although  these  rates  are  rather  low  wo  still  anticipate 
reasonable  gains  since  the  field  induced  lines  undergo 
collisional  narrowing.  This  phenomenon  occurs  if  the 
molecule  undergoes  many  collisions  while  traveling 
over  distances  of  about  one  wavelength  |S|.  As  a 
result  the  lines  are  pressure  broadened  with  width  well 
helmv  the  doppler  width  at  densities  of  I  amapat. 

lixperimental  and  theoretical  data  on  the  collision 
narrowed  lincwidth  arc  available  for  II-.  from  ref. 

(6 1  we  take  the  value  ot  pressure  broadening  coeffi¬ 
cient 

A(X~  1 )  =  2.1  X  10  A'jamagat |  cm  1 


a.v  ->  1 

=  (5  X  10  ")  cm 

where  AX' '.V  is  the  relative  population  inversion. 
Keeping  this  quantity  fixed  the  gain  is  independent 
of  the  density  .V  in  the  pressure  broadened  regime, 
since  the  lincwidth  is  proportional  to.V, 

Higher  gains  can  be  obtained  by  increasing  the 
strength  of  the  applied  field  I'  For  a  statu  meld  elec¬ 
tric  breakdown  will  occur  at  field  strength  larger  than 
the  value  of  105  V/cm  considered  here.  For  the  alter¬ 
nating  electric  field  produced  by  a  giant  laser  pulse 
much  higher  values  arc  possible.  For  example  a  field 
strength  of  I05  V/cm  would  correspond  to  a  laser  in¬ 
tensity  of  2  X  10s  W'cnr .  but  intensities  of  up  to 
10ll)  W/cirr  have  been  achieved  for  picoseconu  pulses 
without  optical  breakdown  [9], 

In  order  to  establish  the  feasibility  of  gain  switching 
using  a  high  power  laser  pulse  as  the  source  ot  out  elec¬ 
tric  field,  we  use  the  following  expression  for  (he  mix¬ 
ing  of  lasing  and  intermediate  states  a  and  r  due  to 
our  switching  laser  signal  |20| 


I  a( /))  =  |a>  + 


<rl/W(t||a> 

-  /:'  +  hvn 


Ii' 


■>‘o' 


(rlfl-fX-'la) 

\  ~  '■  r 


In 


"°f  )|r>. 


(o) 


where  F’*’'  denotes  the  positive  and  negative  narts  of 
the  injected  field,  and  i>n  is  the  frequency  ot  the  tntcct- 
cd  field. 

From  eq.(6)we  see  that  the  inicctcd  l.i  -er  field 
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mixes  states  :n  much  the  same  way  as  tines  an  ortiinarv 
tie  field.  However  as  noted  above  the  fields  are  now 
much  larger  and  tins  would  increase  the  a  am  futures 
obtained  above  by  at  least  two  orders  01  magnitude. 

In  this  ease  earn  trail  res  ot  several  percent  cm  1  could 
be  expected  tor  the  Ns  molecule  also,  l  or  tins  mole¬ 
cule  elticient  ■  rational  pimipma  techniques  such  as 
electron  itnpa.:  I0|  and  discharges  [  11 1  are  available 

From  the  serrucbssicjl  point  ot  view,  application  oi 
a  time  dependent  electric  Held  at  trequencv  :'n  will 
lead  to  an  induced  transition  dipole  moment  oscillat¬ 
ing  at  liequenctes  r  +  i-'(,  and  r  i  ().  where  r  ts  the 
vibrational  transition  frequency Hence  stimulated 
emission  will  be  observed  at  both  sidebands,  in  a  quan¬ 
tum  mechanical  treatment  we  would  have  to  take  into 
account  the  tact  that  the  pain  cross-section  tor  both 
emission  f  requencies  arc  not  exactly  equal  |ld.le.|x|. 

In  the  present  work  the  notion  of  inducing  a  dipole 
transition  matrix  element  via  a  tune-dependent  laser 
held  was  an  outgrowth  ol  our  calculations  involving 
a  dc  electric  held.  Those  calculations  were  based  upon 
the  earlier  work  of  Overhauser  (20| .  Upon  completing 
our  analysis,  it  became  apparent  that  there  is  a  close 
connection  with  the  earlier  work  of  Harris  et  al.,  in 
which  they  utilized  an  incident  laser  to  induce  transi¬ 
tions  from  metastable  atomic  levels  |1S.I'»|.  A  more 
detailed  discussion  of  the  ptesent  work,  connection 
with  the  work  of  Hams  et  al..  and  the  question  of 
population  inversion  in  our  11  •>  system  will  he  given 
elsewhere. 

The  authors  gratefulW  acknowledge  helpful  discus¬ 
sions  with  T.  Avizoms.  k.  Kompa  and  S.  Rabin. 
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rime-delayed  measurement  of  naturally  broadened  transitions  can  lead  to  a  narrowing 
of  the  linewidth.  Moreover,  under  appropriate  conditions,  it  may  result  in  the  appearance 
of  a  dip  at  the  line  center.  An  analysis  of  time-delayed  measurement  thus  provides  a 
theoretical  basis  for  useful  optical  techniques  yielding  high  spectral  resolution.  Such  an 
analysis  is  presented  in  this  work. 

I.  INTRODUCTION  lay  between  the  system  preparation  and  the  obser¬ 

vation  of  emitted  radiation. 

It  has  traditionally  been  one  of  the  main  en-  Knight  and  Coleman4  have  shown  that  a  similar 

deavors  of  spectroscopists  to  develop  measurement  result  may  be  achieved  in  a  system  of  two-level 

techniques  yielding  ever-higher  resolution.  In  opu-  atoms  weakly  driven  by  an  exponentially  decaying 

cal  experiments  the  precision  of  the  measurement  is  laser  pulse.  In  this  system,  although  the  lower  lev  - 

often  limited  by  the  broadening  of  the  linewidth  el  is  stable,  the  exponential  decay  of  the  pulse  am- 

caused  by  the  interaction  of  the  system  under  in-  plitude  has  the  same  effect  on  the  fluorescence 

vestigation  with  its  environment,  such  as  Doppler  spectrum  as  would  the  exponential  decay  of  the 

effect,  collisions,  and  spontaneous  decay.  It  might  lower  level.  Metcalf  and  Phillips'  have  shown  that 

seem  that,  after  the  Doppler-  or  collision-  despite  the  loss  of  signal  associated  with  time- 

broadened  width  has  been  eliminated  using  one  of  delayed  detection,  it  may  still  prove  very  useful  in 

the  many  schemes  introduced  in  the  past  for  that  a  number  of  applications.  For  example,  as  em- 

purpose.1  the  natural  linewidth  remains  the  ulti-  phasiz.ed  in  Ref.  2.  this  technique  would  allow  us 

mate  limit  to  high-resolution  spectroscopy.  to  measure  the  difference  ( ya  —  ft, )  directly,  and 

Recently,  we'  have  proposed  and  analyzed  some  therefore  to  a  much  higher  precision  than  could  be 

spectroscopic  techniques  which  provide  resolution  obtained  from  independent  measurements  of  y0  and 

beyond  the  natural  linewidth.  These  considerations  yh. 

are  based  on  the  fact  '  that,  in  the  transient  regime.  We  note  that  transient  line-narrowing  spectros- 

the  probability  for  induced  transitions  in  a  two-  copy  has  a  number  of  similarities  with  detection 

level  system  interacting  with  a  monochromatic  schemes  developed  earlier  to  achieve  resolution 

electromagnetic  field  is  not  weighted  by  a  Lorentzi-  beyond  the  natural  linewidth  in  Mossbauer.” 1 

an  of  width  yab  siya  +  yb,  but  rather  s)'a  -  )v  level-crossing/- 11  and  Lamb-shift i: ~ 14  expen- 

[ Ya  and  Yb  are  the  decay  rates  of  the  two  levels.  ments.  The  common  feature  of  these  experiments 

Note  that  in  this  paper,  we  call  ya  and  yb  the  am-  is  to  discard  the  part  of  the  radiation  emitted 

plitude  (rather  than  population)  decay  rates.  Thus  shortly  after  the  preparation  of  the  system  and  to 

}'«  and  yb  are  twice  as  large  as  in  the  usual  nota-  collect  only  the  delayed  and  exponentially  weak- 

tion.]  As  a  specific  example  to  demonstrate  such  ened  signal. 

transient  line  narrowing,  we  have  proposed  an  ex-  In  this  paper  we  present  a  fully  quant  uni- 

perimental  setup,  inspired  from  delayed  detection  mechanical  treatment  of  time-delayed  spectroscopy, 

level-crossing  spectroscopy,  that  utilizes  a  time  de-  both  for  weak  and  strong  incident  fields.  In  the 
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discussion  of  Ref.  2  the  atoms  were  assumed  to  he 
driven  by  a  weak  classical  field,  so  that  a  perturba¬ 
tive  treatment  can  be  used.  In  general,  the  strong- 
field  dynamics  is  substantially  different  from  the 
weak-field  one.  Since  the  time-delayed  spectrum 
depends  sensitively  on  the  temporal  behavior  of  the 
system,  one  might  expect  that  the  inclusion  of 
power  broadening  would  lead  to  a  different  time- 
delayed  spectrum.  In  fact,  as  is  shown  in  the 
second  part  of  this  paper,  it  can  lead  to  the  appear¬ 
ance  of  a  transient  dip  at  line  center.  This  din  at 
line  center  may  prove  to  be  a  useful  technique  to 
determine  accurately  the  position  of  the  transition. 
Thus,  time-delayed  fluorescence  measurements 
have  the  capability  of  providing  high  spectral  reso¬ 
lution.  not  only  through  the  line-narrow tng  effect, 
but  possibly  via  strong-field  "transient  dip  spec¬ 
troscopy". 

The  goals  of  the  present  paper  is  threefold. 

First,  we  show  that  the  results  previously  obtained 
semiclassically  are  recovered  exactly  in  a  fully 
quantum-mechanical  treatment.  Second,  we  extend 
the  previous  "weak-field"  calculations  to  arbitrary 
strength  driving  fields.  Third,  we  give  an  intuitive 
physical  picture  of  transient  line  narrowing  and 
transient  dip  spectroscopy,  based  on  a  well-known 
feature  of  the  Rabi  problem. 

The  remainder  of  this  paper  is  organized  in  the 
following  way:  In  Sec.  [f,  we  give  a  fully 
quantum-mechanical  theory  of  transient  line  nar¬ 
rowing.  considering  an  atom  weakly  driven  by  ei¬ 
ther  cw  or  pulsed  excitation.  In  the  case  of  puised 
excitation,  square  pulses,  as  well  as  exponentially- 
decreasing  pulses,  are  considered.  We  show  that 
this  leads  to  exactly  the  same  results  as  the  semi- 
classical  theory.*'4  provided  spontaneous  emission 
directly  between  the  two  states  under  consideration 
is  neglected,  and  the  emitted  photons  from  the  two 
states  are  distinguishable.  In  Sec.  Ill  the  effect  of 
power  broadening  on  the  time-delayed  spectrum  is 
studied.  It  is  shown  that  the  inclusion  of  power 
broadening  leads  to  the  appearance  of  a  transient 
dip.  Finally,  Sec.  IV  is  a  summary  and  discussion. 
Throughout  the  paper,  natural  units  fi=c  =  l  are 
used,  unless  otherwise  stated. 

II.  QUANTUM-MECHANICAL  THEORY  OF 
TRANSIENT  LINE  NARROWING: 

WEAK-FIELD  LIMIT 

A.  Four-level  atom  driven  by  weak  cw  radiation 

We  first  consider  the  system  previously  investi¬ 
gated  by  Meystre.  Scully,  and  Walther.*  This  sys¬ 


tem  consists  of  an  atom  with  two  unstable  levels  ,j 
and  b  and  a  weak  cw  field  dm  me  the  atom  Horn 
the  lower  level  h  to  the  upper  level  a  If  one  in¬ 
cludes  the  lower  levels  to  winch  a  and  b  decav.  this 
may  be  considered  as  a  four-level  atom  see  I  ic  1 
We  prepare  the  atom  in  level  b  and.  as  the  field 
dm  .*s  the  atom  to  level  a,  we  count  the  photons 
emitted  following  the  a  -<  transition,  starting  a 
finite  time  0  after  the  atom  is  prepared  1  lie 
counting  rate  is  measured  as  a  function  of  the  de¬ 
tuning  between  the  laser  frequency  k  ,  and  lac  en¬ 
ergy  separation  AJfl  between  a  and  b 

Under  the  rotating-wave  approximation  and  in 
the  interaction  picture,  the  wave  function  may  be 
written  as 

1  t/x f  1 )  =  a"! f )  a.t.V  -  1  ik'.T  ~ab'I'  b.  \  k  ■ 

-  tv.V  Uk'.k 

-  f  r/kVe-.1'’  .k\  >  .  2  . 

with 

«.'•«)!>--  />..\  k*. ,  >  .  :: 

Here  .V  is  the  number  of  photons  in  the  driving 
field.  a.t.V  -  1  i k,, >  denotes  the  state  in  which  the 
atom  is  in  the  state  a  and  t.V  -  1 1  photons  of  mode 
kn  are  present  in  the  radiation  field  land  similarly 
for  other  eigenstatesi.  and  kn  is  the  wave  vector  of 
the  driving  field.  k|  and  k-  denote  the  modes  ol 
the  photons  emitted  following  the  a  -*<•  and  b 
transitions,  respectively.  In  writing  Eq.  (2.1i  we 
have  assumed  that  the  energy  separations  between 
any  two  levels  are  sufficiently  different  from  one 
another  that  the  photons  k0,  k|.  and  k;  are  distin¬ 
guishable.  We  also  have  assumed  that  spontaneous 
emission  from  a  to  b  can  be  neglected.  Substitut¬ 
ing  Eq.  (2.1)  into  the  Schrodinger  equation,  we  ltn- 
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FIG.  I.  Four-level  atom  driven  bv  cw  radiation. 
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mediately  obtain  the  following  equation-  tor  the 
probability  amplitudes 

i--aat/i  -  v  V  V/ *•  exp( / 1 A,,|r]<r'>'rl 

dt  " 

-  f  i/k.|.V/*-  e.\p(<  iAj,.  A,  /|<C-  i  . 


i  (A  t'  i  Y.W.-  exp(  hAj(, 

dt 


f  JAM'  exp|i  A„ 


s  7  <i  ■  i 


lb 


^nr‘k-  U>- M  k-  e\p(  /  <  AJ( 

A  t  ‘t  ]i  r"'1  t  -  . 

>  *  . 
a.  ,'L 

y«k-U)  =  .Mk-  exp[  n  Aw 

A  ,  7  ].  ."  . 

-  2  3d 

lubject  to  the  initial  condition. 

Fq  2  2  or 

«*(0>=  1  . 

2  4.i 

rtll(0l  —  <r‘k-  lOl-o^  K))-  0 

2  4b 

A,;  is  the  energy  separation  between  ihe  states  / 
and  j.  A,  is  the  energy  of  a  photon,  i.e..  A.  -  A 
and  M k-  is  the  atomic-transition  amplitude  defined 
as 


where  i  a  or  />  and  ;  i  or  J  lor  the  preset!!  ease 
Hi  denotes  the  interaelion  Hamiltonian  in  the  in- 
teraetion  picture 

III  order  to  have  an  unambiguous  definition  .it 
the  tune -delated  power  spectrum,  one  needs  to  in¬ 
troduce  a  model  detec-tor  into  the  problem  In 
this  paper,  we  consider  the  same  scheme  as  m  Ret 
i  e  .  we  detect  the  total  number  01  photons  -pot. 
laneousls  emitted  following  the  transition  a  •<  as 
a  iuiiciioii  oi  (he  detuning  A  between  the  driving 
field  and  the  atomic  energy  separation 
AJfr  A  A  AJ(1  We  then  define  the  tune 
delavcd  spectrum  V  i  A.tC  as 
r  x  , 

N.A.tt'  -YaJ,  J,i  tl“'1  i1  ’•  -ha 

i.e..  as  the  number  of  photoeounls  from  time  ti  on 
This  can  be  reexpressed  as  the  number  of  emitted 
photons  at  t  =  x  minus  that  it  r  =  0. 

.V ( A .())  =  f  dk,  ixl 

-  f  d k,  a{-  U  =6)  :  .  2.6b 

Thus,  the  time-delayed  photocounts  can  be  ob¬ 
tained  by  solving  Eqs.  (2.3)  and  substituting  the 
solutions  into  Eqs.  (2. but  or  (2.6b1,  Details  of  the 
calculation  are  shown  m  Appendix  .A  for  the 
weak-field  limit  Here  we  only  show  the  result 


(j.k,Hi  i,0)  =.Wv-exp[  — ttA,, -A  i/)  .  1 2.5 


2.V  M 

.  s  K0  14 

.V(A.0)  = - , - ; - 


expt-2)'udi  expi  -  2;><K 
~  2r.~  T~  -Y„ 


exp,'  '  I1,0 '  I A  sinAb  -  yab  eosAR ' 

*  }'aS 


where  bab  =y„  - yb ,  yab=ya+Yi,.  and  ya  and  yb 
are  the  decay  rates  of  levels  a  and  b,  respectively. 

We  note  that  in  the  limit  6—*0,  Eq.  (2.7)  yields 
the  usual  Lorentzian  of  width  yab-  However,  if  0 
is  sufficiently  large,  only  the  first  or  the  second 
term  in  the  large  parentheses  of  Eq.  (2.7)  remains, 
depending  upon  which  of  ya  and  yb  is  largest. 

Thus,  the  dependence  of  /V(A,0)  on  A  is  determined 
mainly  by  the  Lorentzian  prefactor,  and  the  line- 
width  approaches  the  difference  bab.  This  is  the 
transient  line-narrowing  effect.  In  fact,  Eq.  (2.7)  is 
identical  with  the  expression  obtained  earlier2  using 
the  density-matrix  equations  and  a  classical 
description  of  the  field.  Here,  Eq.  (2.7)  is  obtained 
using  j  fully  quantum-mechanical  approach.  The 
u-irmalized  power  spectrum  A’(A.f))  is  shown  in 
Fig.  2  for  y„  =  3  and  yb  —  I ,  and  in  Fig.  3  for 
Ya  =  I  01  and  }'*  =  !. 


O.,.,- ......  i 


FIG.  2.  Normalized  photocounts  .Vi  A.(b  as  a  func¬ 
tion  of  A  for  various  values  of  the  delay  R  =  0.  0.67. 

1.33.  and  2.0  for  the  eases  yd  -  3  and  yb  =  1.  0  is  in 
units  of  ;•(,  ’  If  Ihe  system  under  consideration  is  the 
one  shown  in  Fig.  I.  A  is  the  energy  detuning  between 
the  driving  field  and  the  atomic-energy  separation.  If 
the  system  is  that  shown  in  Fig.  4.  A  is  the  energy  de¬ 
tuning  between  the  central  frequency  of  the  driving  pulse 
and  the  atomic-energy  separation,  and  is  the  decay 
constant  of  the  exponentially  decreasing  pulse. 
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FIG.  V  Same  as  Fig.  2  except  that  (>  —  0.1,2, 4  and 
JV  1-01. 

B.  Three-level  atom  driven  by  weak-pulsed  radiation 

Knight  and  Coleman4  have  pointed  out  that 
transient  line  narrowing  can  also  be  realized  in  a 
system  of  two-level  atoms  driven  hv  a  weak  ex¬ 
ponentially  decreasing  pulse.  In  fact,  the  subnatur¬ 
al  linewidth  observed  in  Mbssbauer  experiments" 
involves  basically  the  same  system.  For  conveni¬ 
ence  we  introduce,  in  addition  to  the  two  levels  a 
and  b,  a  third  level  c  to  which  the  upper  level  a 
can  decay  (see  Fig.  4).  Comparison  of  Figs.  1  and 
4  shows  the  strong  analogy  between  this  system 
and  that  discussed  previously.  The  only  difference 
is  that  the  exponentially  decreasing  pumping  rate 
arises  now  from  the  exponential  decay  of  the  pulse 
amplitude  in  Fig.  4  whereas  it  came  from  the  ex¬ 
ponential  decay  of  level  b  in  Fig.  1. 

We  describe  the  radiation  via  its  spectral  ampli¬ 
tude  k  I.  For  an  exponentially  decaying  pulse, 
this  would  be  a  l.orentzian.  However,  we  keep 
<5U)  arbitrary,  so  that  our  result  is  valid  for  any 
pulse  shape.  Let  us  introduce  the  operator  B 
as1"'17 

fl7=  /  rfk<5(k)a,*  ,  (2.81 

where  a  k‘  is  the  usual  photon  creation  operator. 


PULSE 


b 


FIG.  4.  Three-level  atom  driven  by  radiation  pulse. 


The  pulsed  radiation  can  then  be  represented  by 

•V* )  —  —A—- 1  B*)s  vacuum  )  ,  i2.9i 

v  .V ! 

where  .V  is  the  number  of  photons  in  the  pulse. 

For  simplicity,  we  take  .V  =  1  and  consistently  lim¬ 
it  ourselves  to  the  weak-field  limit. 

As  before,  we  write  the  wave  function  as 

;  t/’U )  >  =  i )  |  u,0  >  -t-  f  d  k'o t)  l  b,  k  > 

-+-  f  £/k,acjMii  c.k'j)  .  (2.10i 

with 

tf(0)>  =  J  d  k<4(  k  )  b.k>=  b,l*>.  (2.11) 

Substituting  Eq.  (2.10)  into  the  Schrodinger  equa¬ 
tion,  we  obtain 

i-  «‘,<rl  =  J  dkM  exp[/ ( Aa()  -  A  )r]«k-(t) 


-  f  d  k  (A/  *•  exp[i  (  Aar  -  k  ,  u]nfk*  ( f '  . 

i2.12al 

i  dr 17 v’tf  1  =  -W T e\p[  —  i  ( Aa4  —  k )/]«a(t)  ,  (2. 1 2b> 

i <rkMf  1  =  ,Vf  exp[  —  i  i  Aar  —  k  \  ) t\aa( 1 1  .  <2. 1 2c 1 

where,  as  before,  we  have  assumed  that  \ab  and 
Aar  are  sufficiently  different  from  each  other  so 
that  the  photons  k  and  k|  are  distinguishable. 
Equations  (2.12'  are  subject  to  the  initial  condition. 
Eq.  (2.1 1),  or 

<r  ki0)  =  <5(  k  1.  aa(0)  =  «k*  (0)  =  0  .  (2.13' 

The  number  of  photons  emitted  following  the 
a-*c  transition  from  the  time  t  =0  on  is  given  by 

iV(A.(ri=2;-a  f*dt]  rta(f,'  ; 

=  J"d k  i  ack‘ it  =  x  ); '  —  j d k* i  ti kMt  =  0 '  *  . 

(2.14) 

This  can  be  obtained  by  solving  Eqs.  (2.12)  and 
substituting  the  solution  into  Eq.  (2. 14'.  Details  of 
the  calculation  are  shown  in  Appendix  B.  neglect¬ 
ing  spontaneous  decay  from  a  to  b  t;  a  is  the  decay 
rate  of  the  a  ~i  transition).  .Y(A.(Jl,  of  course, 
depends  on  the  pulse  shape  because  the  solutions  of 
Eqs.  (2.12)  are  different  for  a  different  choice  of 
<J(  k ). 

For  the  case  of  an  exponentially  decaying  pulse. 
.YiA.O)  is  again  given  by  Eq.  (2.7),  Therefore. 
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Figs,  2  and  3  also  give  the  normalized  count 
.Y(A,tM  lor  the  system  considered  here.  We  note, 
however,  that  kn  is  now  to  be  interpreted  as  the 
central  frequency  of  the  pulse,  A  =  k„  —  Aai),  and  yh 


is  the  decav  constant  of  the  pulse  rather  than  of 
the  atomic  level. 

For  the  case  of  a  square  pulse  of  duration  we 
obtain 


,V(A,0>  = 


~  :  8s  Ya 
A2  +  Ya 


in  —  6  — 


expt  -2yar, ,  I  -  ex  pi  -2yat» 

-Ya 


2  expt  -  ?'„/,)  H  ya  cosA  f„  -  A  sinAf,,  i  -  2  expt  -ya0^ya  eosAO  -  AsinAfU 

A:  -k  v,: 


i  g,  -  expl  -  2ya0) 
±  r  Ya 


[  1  4-  expt  2)Vu )  -  2  exptya/0 )  cosA/,,] 


:  i2.15.i 

for  0.  i  . 
for  0  ■s  . 
<2. 1 5b> 


g,  is  a  constant  depending  on  the  pulse  parameters,  k0  is  again  the  central  frequency  of  the  pulse,  and 
A  =  k„  —  ±ab.  The  normalized  count  ,Y(A,(>)  is  shown  in  Fig.  5  for  t„=  1  and  ya  =  3,  and  in  Fig.  b  for  t„=  1 
and  ya  =  1.01.  ,V(  A.fl)  takes  the  same  form  for  all  0\  exceeding  /„•  This  is  because,  for  ■>  t„,  the  situation 
is  analogous  to  making  time-delayed  measurement  on  a  two-level  system  with  the  levels  a  and  c.  We  see 
that  the  line  narrowing  does  occur  with  a  square  pulse  although  the  narrowing  is  not  as  dramatic  as  that  as¬ 
sociated  with  an  exponentially  decreasing  pulse. 


C.  Physical  interpretation  of  transient  line  narrowing 

Since  transient  line  n arrowing  arises  from  the 
behavior  of  the  system  in  the  transient  regime,  it  is 
natural  to  study  the  temporal  behavior  of  the  sys¬ 
tem  in  order  to  better  understand  and  interpret  the 
effect.  Intuitively,  it  is  not  difficult  to  understand 
the  line-narrowing  effect.  It  is  based  on  the  well- 
known  (but  perhaps  anti-intuitivei  fact  that  the 
Rabi  frequency  is  larger  for  larger  detunings,  and 
that  therefore  the  excitation  and  the  depletion  of 
the  population  of  the  upper  level  (level  a  in  our 
system)  are  slower  when  excited  exactly  on  reso- 


% 

V\ 
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FIG.  5.  Normalized  photocounis  .ViA.t))  as  a  func¬ 
tion  of  energy  detuning  for  various  values  of  the  delay 
0  0.  0.5.  and  ()■>  I  for  the  case  y„  =3.  The  driving 
field  is  assumed  to  be  in  the  form  of  a  square  pulse  of 
duration  r„  -  I.  0  is  in  units  off,,. 


nance.  This  means  that  the  remaining  population 
after  some  delay  time  0  is  relatively  a  large  number 
when  excited  on  resonance,  thus  leading  to  the  line 
narrowing.  In  this  section  we  show  quantitatively 
that  the  above  interpretation  is  indeed  true. 
Although  our  discussion  here  is  limited  to  the  case 
of  an  atom  driven  by  pulsed  radiation,  it  can  equal¬ 
ly  be  applied  to  the  system  considered  in  Sec.  II  A 
provided  that  the  appropriate  redefinition  of  k„. 
yb.  etc.,  is  made. 

The  integrand  n^U)  *  in  Eq.  (2.14)  for  the  pho¬ 
tocount  is  the  probability  for  the  atom  to  be  in  the 
upper  state  a  at  time  r.  For  the  case  of  an  ex¬ 
ponentially  decreasing  pulse,  we  find  using  Eq. 

(B4)  of  Appendix  B, 


FIG.  t>.  Same  as  Fig.  5  except  that  1.01 
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FIG.  7.  Probability  P"  (A./ 1  as  a  function  of  time  r  for 
two  values  of  energy  detuning  A  -  0  and  2  for  the  cases 
=  3  and  >'*  =  I  (exponentially  decreasing  pulse).  The 
time  is  ii.  Jims  of  )'»  l. 


pa(\,t)~  j  aa(t)  j  ’ 

I  Se  !  ^ 

=  ~ ( exp<  -  2ya r  I  +  exp(  -  2yb 1 1 
A-  +6a*(, 

-2expl  —  )'otf)cosAr]  .  (2.16) 

p“(A,/)  is  plotted  in  Fig.  7  as  a  function  of  time 
for  two  different  values  of  detuning  A  =  0  and  2  for 
the  ease  y„  =  3  and  yb  =  1 ,  and  in  Fig.  8  for 
ya  =  1.01  and  yb  =  1.  pa(A,t)  initially  increases  as 
the  driving  field  pumps  the  atom  to  level  a,  goes 
through  a  peak  and  eventually  decreases  to  zero  be¬ 
cause  of  spontaneous  emission.  The  peak  occurs 
earlier  and  has  a  smaller  value  for  A  =  2.  This  re¬ 


FIG.  8.  Same  as  Fig.  7  except  that  ya  =  1.01. 


fleets  the  above-mentioned  fact  that  the  Rahi  fre¬ 
quency  is  larger  for  larger  detumne'  nd  is  ihc  key 
point  in  understanding  the  line-narrowing  effect. 
According  to  Eq.  (2.14),  A(A,0)  is  the  area  under  the 
probability  curve  between  r  =  6  and  x .  We  see 
from  Fig.  7  that,  as  6  is  increased,  this  area  de¬ 
creases  faster  for  larger  A.  This  is  even  more  clear 
in  Fig.  8.  The  integral  of  pa(A,t)  between,  for  ex¬ 
ample.  I  ~2  and  x  is  only  a  small  fraction  of  the 
area  between  t  -0  and  *  for  the  case  A  =  2.  while 
it  is  still  a  large  fraction  for  the  case  A  =0.  Thus, 
the  larger  A.  the  faster  the  system  emits  spontane¬ 
ously  and  leaves  a  smaller  number  of  photons  to  be 
emitted  after  some  delay  time  0.  This  directly 
leads  to  the  narrowing  of  the  hnewidth. 

For  a  square  pulse  of  duration  the  probability 
takes  the  form 


p“t  A.l)  — 


\  g  ;  - 

i  a“U)  |  :=  — ~ — 1  +  exp(  —  2yat)- 2exp( -}-8r)cosAt],  for  t  <t0  (2.17a) 

+  y; 

,  |g,  |Jexp(— 2}'„/) 

|  ot®( f )  |  ■  = - — - , - [exp(2:Jt„)-i-l  -  2 exp( }•„/,, )cosA/0],  for  im(,  (2.17b) 

*  +VJ 


where  Eqs.  (B8)  of  Appendix  B  have  been  used.  This  probability  is  plotted  in  Fig.  9  as  a  function  of  time 
for  two  different  values  of  detuning  for  the  case  tl(  =  I.  ya  -  1.01.  Here  again  we  see  that  the  peak  occurs 
earlier  and  has  a  smaller  value  for  larger  detunings,  although  the  effect  is  not  as  strong  as  before.  This  ex¬ 
plains  a  relatively  weak  li-  •  narrowing  with  a  square  pulse. 


III.  STRONG-SIGNAL  REGIME:  TRANSIENT  DIP  SPECTROSCOPY 

Up  to  now,  we  have  restricted  ourselves  to  the  weak-field  limit,  in  which  the  "««*  of  perturbative  treat¬ 
ment  is  justified.  We  now  relax  this  restriction  and  study  time-delayed  spectroscop.  ,i  the  strong-signal  re¬ 
gime. 


i 
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The  analysis  departs  from  that  of  the  weak-field  limit  in  that  instead  of  using  Eq.  (AS)  of  Appendix  A,  we 
now  solve  Eqs.  (A2)  and  (A3)  exactly  for  a“U )  and  ab(t).  These  equations  yield  a  second-order  differential 
equation  for  aaU)  [and  also  for  a*(f)J,  with  the  solution3 

a#U)  =  —  tv  N M exptu.t)  JQ  dt\  exp[(u  |  —  tt; )/,  ]  , 

(3.1) 

where 

U,  =  l  —Yab  —  »A  +  [(6a*  +  /A)2  —  4A  |  M ^  - ] 1/2 1  /2 

,  (3.2a) 

"2  =  1  -Fa* -1  A  — [(&„*+, A)2-4A  |Afroj2],/2|/2 

(3.2b) 

Using  Eqs.  (3.1)  and  (2.3c)  we  get 
,  _  r< 

Jg  dtK  exp[u2  —  i  ( Aaf  —  k  t  ]r ,  Jg  dt,  exp(u  i  —  u,  )l;  ■ 


(3.3) 


The  number  of  photons  emitted  during  the  time  interval  (0, oc  )  following  the  a—*c  transition  is  again  given 
by  Eq.  (2.6a)  or  (2.6b).  After  tedious  but  straightforward  algebra,  we  obtain 

<V(A,0)  =  [4 yaN  |  |  2exp(  -y<*fl)]/p 

yab  cosh[ 0 Vp cos( <f>/2 ) ]  +  v'p cos( 6/2 1  sinh[  0  v  p cos( 6/2 ) ] 


Yab  — pcos:(<A/2) 

Yab  cos[ 0 v  p  sin(<£/2 )]  —  v'p  sin! 6/2 )  sin[ 0Vp  sin( 6/2 ) } 


Yab+P  stn-(<6/2 ) 

where 

p  =  [(6;*-A:-4A  IMrjJ)J+4«iA:]l/J  . 
and  6  is  determined  by 

cos<*  =  (6j*  -A2-4/V  |  M j,-  |  2)/p  , 

sin<^  =  26a*A/p  . 


(3.4) 

(3.5' 

(3.6a) 

(3.6b) 


Equation  (3.4)  is  the  main  result  of  this  section.  The  time-delayed  count  A(A.0)  is,  in  general,  a  complicat¬ 
ed  function  of  the  system  parameters,  and  we  first  consider  some  limiting  cases. 

f - 


A.  Weak-field  limit 

In  the  weak-field  limit  ,V— »0,  we  obtain 
p-6;h  +A:  , 

,r-  ...... - 


(3.7a) 


b  / 


cos<£~(6;*  -  A :>/p  ,  (3.7b) 

sin<J~26a*A/p  .  (3.7c) 

Equation  (3.4)  reduces  then  to  Eq.  (2.7)  of  Sec. 

II  A,  as  it  should. 

B.  No  time  delay 

In  the  limit  0*0  (no  time  delay),  Eq.  (3.4)  be¬ 
comes 

•v  U/n  ::r*b 

•V(  A,0  -0)  -  -  -;  -7 - r  —  — 7 — - - 

)  /»'  ^  't-  )  ah  f  i,*  Yab  ^ Yu  Yb  ' 


FIG.  9.  Probability  P“  as  a  function  of  time  /  for  two 
values  of  energy  detuning  A  0  and  2  for  the  case  r„  -  I 
(square  pulsei  and  y,  -  1.01.  /  is  in  units  of  f„. 


1 3.8' 


which  is  the  well-known  power- broadened 
Lorenman  line  shape. 
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C.  General  case 

We  now  return  to  the  general  formula  Eq.  (3.4). 

In  Figs.  10  and  11,  Y(A,0)  is  plotted  as  a  function 
of  the  detuning  A {=kn-&ab )  for  different  delay 
times  8  and  for  different  field  intensities  for  the 
case  ya  =3  and  yb  =  1.  We  immediately  see  a  new 
feature  of  the  power-broadened  spectrum,  namely, 
the  appearance  of  a  dip  at  line  center  for  large  8 
and  for  high  enough  field  intensities.  A  close  in¬ 
vestigation  indicates  that,  for  a  fixed  value  of  the 
field  intensity,  and  as  0  is  increased  from  zero,  the 
linewidth  first  decreases  (transient  line  narrowing). 
However,  as  8  is  further  increased,  the  line 
broadens  back  until  the  linewidth  becomes  roughly 
the  same  as  that  at  0  =  0.  If  8  is  further  increased 
beyond  this  point,  a  narrow  dip  appears  at  line 
center,  and  becomes  deeper  and  wider  as  8  is  fur¬ 
ther  increased.  Since  it  always  appears  at  line 
center,  it  may  help  locate  the  center  of  a  transition 
with  improved  accuracy  over  standard  methods. 

The  dip  appears  with  a  smaller  delay  for  stronger 
fields.  Thus,  one  can  operate  without  significant 
loss  of  signal.  For  example,  for  a  Rabi  frequency 
G(  =  Y  i  Mi'  |  *)=  1.5,  the  dip  appears  already  at 

8  =  2.  On  the  other  hand,  for  G  =0.6,  the  dip  still 
does  not  appear  even  for  8  as  large  as  4.  In  Fig. 

12  we  show  iV(A,0)  at  a  fixed  6  as  G  is  varied. 

Here  again,  we  see  clearly  that  the  dip  appears 
more  easily  for  stronger  fields. 

In  Fig.  13  we  show  the  same  ;V(A,0)  versus  A 
curve  for  the  cases  ya  =  1.01  and  yb  =  1,  and  G  =  1 
for  various  8.  We  immediately  notice  the  strong 
oscillatory  behavior  exhibited  in  this  case. 

A  useful  parameter  to  determine  the  behavior  of 
.V(A,0)  is  the  ratio  of  the  Rabi  frequency  to  bab,  or 

r  =  G-/blb  .  (3.9) 


-10  ’•>  '  ’  ”  6“  ’  ~  5  i0 

Ewqy  detuning  A 


FIG.  10.  Normalized  photocounls  Y(A,0l  as  a  func¬ 
tion  of  enerev  detuning  A  for  various  values  of  the  delay 
0  =  0,2,4,6,8,10  for  the  case  ya  =3,  yb  =  I,  and  the  Rabi 

frequency  G  =  Y  t  M  r*  2  -0. 6.  0  is  in  units  of  yZ '. 
o 


i-i  s 


tnerqy  detqn.nq  A 

FIG.  11.  Same  as  Fig.  10  except  that  0  =  0.1,2.4,10 
and  G  =  1.5. 

This  ratio  determines  whether  the  arguments  of  the 
hyperbolic  functions  are  larger  or  smaller  than 
those  of  the  sinusoidal  functions  in  Eq.  (3.4),  and 
therefore  determines  the  characteristic  behavior  of 
Y(A,B).  For  the  case  of  Fig.  13,  r  is  very  large 
(r  =  104)  and  the  contribution  from  the  sinusoidal 
terms  is  important.  As  a  result  Y  ( A,0)  exhibits  a 
strong  oscillatory  behavior.  On  the  other  hand,  r 
is  relatively  small  (r  <  1 )  for  the  cases  shown  in 
Figs.  10  and  11,  and  no  strong  oscillatory  behavior 
is  exhibited  there. 

In  a  sense  the  transient  dip  is  an  opposite 
phenomenon  to  the  line  narrowing  because,  in  or¬ 
der  for  the  dip  to  appear,  the  remaining  population 
of  the  upper  level  a  after  the  delay  time  8  should 
be  a  relatively  small  number  when  excited  on  reso¬ 
nance.  That  is,  the  depletion  of  the  population  of 
level  a  should  be  faster  when  excited  on  resonance. 
That  this  is  indeed  the  case  when  8  is  sufficiently 
large  can  be  seen  in  Fig.  14,  where  the  probability 
paiA.l)=  |  a“{t)  J :  for  the  atoms  to  be  in  level  a  is 
plotted  as  a  function  of  time  for  two  different 
values  of  detuning  for  the  cases  ya  =3.  yb  =  1,  and 
G  =0.6.  As  we  have  already  noted  in  the  previous 
section,  the  peak  occurs  earlier  and  has  a  smaller 
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FIG.  12.  Normalized  photocounts  .VIA .0)  as  a  func¬ 
tion  of  energy  detuning  A  for  various  values  of  (he  field 
intensity  G  =0.3,  0.6,  0.9,  and  1.5  for  the  cases  y,  =3, 
yb  =  I,  and  for  a  fixed  delay  time  (1  =  4.  8  is  in  units  of 
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FIG.  13.  Normalized  photocounts  .V(A,0)  as  a  func¬ 
tion  of  energy  detuning  A  for  various  values  of  the  delay 
tf  =  0,2,6  for  the  cases  ya  =  1.01.  yb=  1,  and  G  =  1.  6  is 
in  units  of  yb'. 

value  for  A  =  2.  This  of  course  accounts  for  the 
narrowing  of  the  linewidth  for  small  values  of  the 
delay  time  0.  However,  when  0  is  sufficiently 
large,  the  A  =  2  curve  final)'-  catches  up  with  the 
A=0  curve.  This  behavior  can  be  understood  if  we 
note  that  for  a  hrge  detuning  the  atom  is  ineffi¬ 
ciently  pumped  t;  level  a  and  therefore  still  a  large 
number  of  atoms  re  avai  able  to  be  pumped  at  a 
large  time.  This  •leans  a  slower  depletion  of  level 
a  at  sufficiently  .;  te  '.tine  when  excited  off  reso¬ 
nance  and  icsds  aircxtly  to  the  appearance  of  a  dip 
at  line  center. 

IV.  SUMMARY  AND  DISCUSSION 

Under  appropriate  conditions,  time-delayed  ob¬ 
servation  of  the  radiation  emitted  during  an  atomic- 
transition  can  lead  to  the  narrowing  of  the 


Time  t 

FIG  14.  Probability  P “  as  a  function  of  time  f  for 
two  values  of  energy  detuning  A  - 0  and  2.  for  the  case 
y,  =  3,  Yh  -  1.  and  G  =0.6.  /  is  in  units  of  yh  For 
is 4,  In/10  is  plotted  as  a  function  of  time  /  and  the  cor¬ 
responding  scale  is  shown  on  the  right  side  of  the  graph. 


linewidth  or  the  appearance  of  a  dip  at  line  center. 
The  linewidth  is  limited  in  the  small-signal  regime 
by  the  difference  in  the  transition  rates  involved. 
This  provides  an  optical  technique  that  yields  high 
spectral  resolution. 

A  dip  may  appear  at  line  center  as  a  result  of 
the  transient  behavior  of  a  system  subject  to  power 
broadening.  Unlike  the  line-narrowing  effect,  the 
width  of  the  dip  is  not  limited  by  (y0  —  yj, ),  and/or 
the  delay  6.  It  can  be  as  small  as  one  likes  if  the 
intensity  of  the  driving  field  is  close  to  the  thresh¬ 
old  for  its  appearance.  This  dip  suggests  a  means 
of  enhancing  resolution  involving  homogeneously 
broadened  systems  while  the  Lamb  dip  is  available 
to  study  inhomogeneously  broadened  systems. 

The  narrowing  of  the  linewidth  or  the  appear¬ 
ance  of  the  dip  is  achieved  at  the  expense  of  some 
loss  of  the  signal.  A  recent  analysis5  indicates  that, 
despite  the  signal  loss,  time-resolved  line  narrowing 
is  highly  desirable  in  a  large  number  of  cases.  If 
one  tries  to  locate  line  center  in  the  absence  of 
spectral  complications,  i.e.,  if  we  know  we  have 
only  one  line  in  the  region  of  interest,  then  it  may 
be  best  to  let  the  delay  time  0— »0  and  collect  the 
maximum  number  of  counts,  since  the  determina¬ 
tion  of  the  line  center  improves  as  the  square  root 
of  the  intensity  of  the  signal.  However,  if  there  are 
complicating  circumstances,  such  as  overlapping 
lines,  etc.,  it  may  be  better  to  use  finite  delay  times 
with  the  attendant  line  narrowing  and/or  dip. 

We  finally  note  that  the  results  obtained  here  via 
a  fully  quantum-mechanical  theory  take  exactly  the 
same  form  as  those  previously  derived  semiclassi- 
cally.  The  origin  of  this  exact  agreement  lies  in 
the  neglect  of  spontaneous  decay  from  level  |  a  >  to 

|f>>. 
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APPENDIX  A:  DERIVATION  OF  EQ.  (2.7). 

THE  CASE  OF  cw  RADIATION 

In  order  to  derive  Eq.  (2.7),  we  first  solve  Eqs. 
(2.3),  From  Eq.  (2.3c)  we  have 
rt 

«kV'>=  -'Me,  J„di i  «P[  -» < A*.  -*i  >f|]«#U,  i  • 

(AD 

Substituting  Eq.  (Al)  into  Eq.  (2.3a).  we  obtain 
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aaU'  =  -iv  A'A/r  exp( -i\t)abU)  . 

(A2) 

Similarly,  from  Eqs.  (2.3d)  and  (2.3b),  we  obtain 
«*(/'=  — iVJVM.-  exp(/A/)aa(f )  . 

lu 

(A3) 

The  decay  rates  ya  and  yb  in  Eqs.  (A2)  and  (A3) 
are  defined  as 

rfl=  |A/r|-’]k=v  ,  (A4) 

- J 

at.)(,)=s(_/)VTvMiMri  /odt,  exp[  —  ya  —  i(A0f 


Substituting  Eq.  (A6)  into  (2.6a),  or  Eq.  (A7)  into 
(2.6b),  we  immediately  obtain  Eq.  (2.7)  for  the 
time-delayed  photocounts. 


APPENDIX  B:  DERIVATION  OF  EQS.  (2.7)  AND  (2.15). 
THE  CASE  OF  PULSED  RADIATION 


In  order  to  derive  Eqs.  (2.7)  and  (2.15)  for  the 
case  of  pulsed  radiation,  we  first  solve  Eqs.  (2.12). 
From  Eqs.  (2.12b)  and  (2.12c),  we  have 

a\{t)  =  i(k)  —  iM  k*  fg  dt,  exp[  —  /  (A„a  —  k)t,  ]aa(/| )  , 

(Bll 

ri 

a\-U)=  Jndt  i  exp[  —  /  ( Aac  —  A ,  )r  ( ]cra( r , )  . 

(B2) 


Substituting  Eqs.  (Bl)  and  (B2)  into  Eq.  (2.12a),  we 
obtain 


di+Ya+Ya 


a“(t) 


— —i  J  dkM{-exp[i(Aai,~k)t]6{\i.)  . 

(B3a) 

=  — igU )  ,  (B3b) 


where  ya  and  y'a  are  the  spontaneous  decay  rates 
from  level  a  to  c  and  b,  respectively,  and  Eq.  lB3b) 
defines  the  pulse-shape  function  g{/). IM7  From 
here  on  we  assume  that  ya  is  negligibly  small. 


)'»  is  given  by  the  same  expression  except  that  it  is 
evaluated  at  k  =  AM  instead  of  k  =  Aoc . 

In  the  weak-field  limit  (.V  small,  or  Rabi  fre¬ 
quency  <7  <<  ya,  yb  ).  the  right  side  of  Eq.  (A3) 
may  be  neglected.  This  means  that  the  major 
source  of  depleting  level  6  is  not  the  driving  field 
but  spontaneous  decay  to  level  d.  We  then  have 

afc(t)2:exp(  ~ybt)  .  (A5i 

Substituting  Eq.  (A5)  into  (A2)  we  then  obtain 
aa(r) 

-  - ' v  V  M  *•„  exPl  -  Ya '  I  fn  dt ,  exp[  (  —  bah  -  i  A ), ,  ] 

(At>) 

and  therefore,  using  Eq.  (Al), 


-k,)]t|  J*o‘ A,exp[(-8at-iA)/:]  . 


(A7) 


I - 

From  Eq.  (B3b)  we  then  have 

-i 

oa(t)=  —i  exp(  —yat)  Jg  dr,  exp(yar,  )g(r, )  . 

iB4! 


and  from  Eq.  (B2)  we  have 

« V,  1/  >  =  (  -  n:M  j  dt ,  expf  -  ya  - 1  ( Aac  -  k ,  >)f , 

r'i 

X  J  dt:exp(yat:  tgIM  . 


The  time-delayed  photocounts  A’  ( A,0)  can  be  ob¬ 
tained  by  substituting  Eq.  (B4)  or  (B5)  into  Eq. 
(2.14).  .V  (A,0)  depends  on  the  pulse  shape  through 
the  function  glr). 

For  an  exponentially  decreasing  pulse,  the  pulse- 
shape  function  g(t)  takes  the  form 

g(f)=gr exp[( -yb~i A)/]  ,  (Bb) 

where  g,  is  a  constant,  A  =  A„-Aa*,  A„  is  the  cen¬ 
tral  frequency  of  the  pulse,  and  yb  is  the  decay 
constant  of  the  pulse.  Substituting  Eq.  iBb)  into 
Eq.  (B4)  or  (B5).  and  using  Eq.  (2.14),  we  obtain 
Eq.  (2.7)  for  the  time-delayed  count  .Y(A,d). 

For  a  square  pulse  of  duration  f„,  we  have 


gs  expt  —  iAi  ),  Ovtv.t„ 

SV"=  |0. 

,.  K0r,  ,,  a  constant,  An  is  again  the  central  fre- 
queii.^  ol  the  pulse,  and  A  =  A„  —  Aab.  Substituting 
Eq.  (B7)  into  Eqs.  (B4)  and  (B5).  we  obtain 


c 1  ct[y>  I 

+  J,  rfliexpl-ra-itA*.-*,)]^  J  dt2c\p[(ya  —/A)/;]  1  B4i 

•*  *o  U  1 

Substituting  Eq.  (B8)  or  (B9)  into  Eq.  t2.14),  we  obtain  Eqs.  (2.15)  for  the  time-delayed  photocounts. 
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The  basic  laws  of  the  phenomenological  theory  of  radiative  energy  transfer  are  derived,  under  certain  conditions,  with¬ 
in  the  tramew  ork  ol  the  stochastic  scalar  wave  theory.  An  equation  of  radiative  energy  transter  is  derived  tor  wave  propaga¬ 
tion  in  a  statistically  quasihomogeneous  medium.  Out  results  relate  the  extinction  and  scattering  coefficients  [which  are 
introduced  heunstically  in  the  conventional  theory  of  radiative  energy  transfer)  to  the  stochastic  characteristics  oi  the 
medium. 


1.  Introduction 

The  subject  of  radiative  energy  transfer  through 
stellar  atmospheres  and  through  turbulent  media  is 
generally  treated  on  the  basis  of  a  phenomenological 
theory  [  1 ,2] .  In  recent  years,  many  attempts  have 
been  made  towards  providing  a  satisfactory  basis  for 
the  conventional  theory  and  to  deliminate  its  range 
of  validity  [3-14] .  These  attempts  have  met  so  far 
with  a  limited  success. 

The  central  quantity  of  the  phenomenological 
theory  is  the  so-called  specific  intensity  of  radiation 
I(R.s)  which  satisfies  the  radiative  energy  transfer 
equation  of  the  form 

s-VrI(R,s)  =  -a(R,  s)I(R.s) 

+  Jd(R,s,s')l(R,s')dn^  +  D(R,s).  (1.1) 

In  this  equation  the  left-hand  side  represents  the  rate 
of  change  of  the  specific  intensity  along  the  s-direc- 
tion.  The  first  term  on  the  right  represents  "extinc¬ 
tion"  due  to  absorption  and  scattering,  the  second 
term  represents  a  contribution  due  to  scattering  from 
all  directions  (the  integration  extending  over  the  com¬ 
plete  4n -solid  angle  generated  by  the  unit  vectors') 
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and  the  last  term  represents  the  effect  of  sources.  The 
functions  a(R,s),  p(R.s,s  )  and  D(R.s)  are  the  so- 
called  extinction  coefficient,  the  differential  scatter¬ 
ing  coefficient  and  the  source  function  respectively. 

In  the  usual  heuristic  model,  the  specific  intensity 
/(R,s)  is  treated  as  a  radiometric  quantity  and  the 
equation  of  radiative  energy  transfer  is  derived  by  an 
intuitive  quasi-geometrical  argument  involving  balance 
of  radiant  energy,  without  elucidating  the  microscopic 
meaning  of  the  extinction  and  scattering  coefficients 
and  of  the  source  function. 

In  this  paper  we  investigate  the  foundations  of  the 
theory  of  radiative  energy  transfer  and  derive  an  equa¬ 
tion  of  radiative  energy  transfer  in  a  randomly  fluc¬ 
tuating  medium  on  the  basis  of  scalar  wave  theory. 

Our  results  elucidate  the  relationship  between  the 
extinction  and  scattering  coefficients  and  the  stochas¬ 
tic  characteristic  of  the  medium. 


2.  Relationship  between  the  specific  intensity  and 
the  coherence  function  of  the  wavefield 

We  consider  a  random  scalar  wavefield  r(r|.  The 
second-order- coherence  function  of  the  wavefield 
is  given  by 

,,2)  =  <p(r|)i'*(r,»).  C.l) 

In  eq.  (2.2),  the  sharp  brackets  denote  the  ensemble 


315 


Volume  37.  number  5 


01»IK  S(  ()\1\U  NR  A  I  IONS 


1  June  i I 


average.  The  coherence  function  l'(r , , r •, ),  when  ex¬ 
pressed  in  terms  of  the  new  variables 


will  be  denoted  by  l'(R,r);  i.e., 

r(fl.r)  (2.3) 

The  ensemble  average  of  the  energy  density  HR) 
and  the  ensemble  average  of  the  energy  tlux  F(R )  are 
related  to  the  coherence  function  of  the  wavefield 
by  the  following  equations  [7| : 

l\R)  =  (1  V)n/?,0)?  f(«)  =  (  1  ik())[7r  l'f/f  .r)|r=(J, 

(2.4M2.S) 

where  c  denotes  the  speed  of  light  in  vacuo  and 
is  the  wavenumber  in  free  space. 

Not  long  ago,  Wolf  [12|  formulated  a  rigorous 
theory  of  radiative  transfer  in  free  electromagnetic 
field.  In  liis  theory,  two  quantities  are  introduced 
which  relate  to  energy  transport:  the  so-called  “an¬ 
gular  components”  of  the  average  electromagnetic 
energy  density  and  of  the  average  Poynting  vector. 

In  analogy  with  these  quantities,  we  introduce  the 
angular  component  of  the  average  energy  density 
Q{R,s )  and  the  angular  component  of  the  average 
energy  flux  T(R,s )  by  the  following  relations: 

U(R)  =  /  Q(R,s)dQ,  F(R)  =  Jt(R.s)  d£2. 

(2.6M2.7) 

The  integrations  in  eqs.  (2.6)  and  (2.7)  extend 
over  the  whole  47r-solid  angle  generated  by  the  unit 
vectors.  If  we  define  the  spectral  density  f(R.K)  of 
the  function  l'(R,r)  with  respect  to  the  variable  r 
by  the  formula 


tive  transfer  for  the  “space  density"  of  radiation 
l  (R)  and  for  the  net  flux  FIR ),  ai  M.  in  iertns  of  the 
specific  !M><*!’.:-i,v  of  radiation  HR.s).  viz. 

i'(R)  =  ^  fllR.s) dll.  (2.111 

F(R)=Js!(R.s)  dll.  (2.12’ 

A  comparison  of  eq.  (2.1  1 )  with  eq.  (2.6)  suggests 
that  the  specific  intensity  of  radiation  HR.s)  can  be 
related  to  the  angular  component  of  the  energy  den¬ 
sity  tJ(R.s)  by  the  -imple  equation 

I(R,s)  =  cQ(R,s).  (2.13 

In  view  of  eq.  (2.9).  we  can  write 


I(R.s)  =  f  f(R.K\ K~  dK.  (2.14) 

0 

The  expression  (2.8)  for  the  spectral  density  j(R. 
JO.  combined  with  eq.  (2.)  4).  gives  a  relationsltip  be¬ 
tween  the  specific  intensity  of  radiation  and  the  co¬ 
herence  function  of  the  wavefield  which  seems  to 
have  been  considered  first  by  Ovchinnikov  and 
Tatarskii  [15] .  Recently,  on  the  basis  of  this  defini¬ 
tion,  rigorous  equations  have  been  derived  by  the 
present  author  for  radiative  transfer  of  energy  and 
momentum  in  free  space  in  the  presence  of  random 
source  distribution  [14] . 

It  is  clear  from  eqs.  (2.7)  and  (2.10)  that,  with 
this  definition  of  the  specific  intensity,  the  usual  flux 
relation  (2.7)  is.  in  general,  not  obeyed.  The  condi¬ 
tion  under  which  the  radiometric  and  field  theoretic 
definitions  of  energy  flux  coincide,  is  given  by 


f{R,K)  =  (2nyi  f  r(R,r)exp{-\K-r)d3K,  (2.8)  J  T(R.s) dll  =  J sl(R.s) dfi>. 


then,  in  view  of  eqs.  (2.4) — (2.7).  we  obtain 


C(R.s)=l  f  f(R,K)  K2  dK. 

c  J 
0 

T(R.s)  =  -  /  f(R,fOK3  dK. 


Equations  (2.6)  and  (2.7)  are  analogous  to  the  i 
pressions  of  the  phenomenological  theory  of  radia- 


On  substituting  from  eqs.  (2.10)  and  (2.14)  for  T  and 
/,  eq.  (2.15)  becomes 

fs(k  k,  )/</?.  AT  ft3*.’ =  0.  (2  16) 

where,  as  before,  s  represents  the  unit  vector  along  A'. 
Consider  a  solution  of  eq.  (2.16)  of  the  form 

/(/},*)  =  /(A.s)fi(A  k());k3.  (2.D 

On  taking  the  Fourier  inverse  of  eq.  (2.8)  and  sub- 
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stunting  for  JiR.K)  from  eq.  (2.17).  it  follows  that 

\)R.r)  =  f  HR  ,s)cxp(iA:0s-r)d<2.  (2. IS) 

This  representation  of  coherence  function  has  been 
adopted  by  many  authors  in  their  investigations  on 
the  foundation  of  the  theory  of  radiative  energy  trans¬ 
fer. 

In  a  recent  paper,  Collett,  Foley  and  Wolf  [lb] 
have  shown  that  if.  for  propagation  in  free  space,  the 
coherence  function  T(R.r)  admits  the  representation 
eq.  (2.18)  then  the  wavefteld  hr)  is  necessarily  statis¬ 
tically  homogeneous.  Furthermore  one  can  establish 
the  following  theorem:  Theorem  2.1:  The  coherence 
function  l'(/f,r)  of  the  wavetield  in  a  statistically 
homogeneous  medium  is  independent  of  R.  (i.e.,  is 
a  function  of  the  difference  variable  r  only)  in  the 
special  case  when  the  incident  field,  which  satisfies 
the  free  space  wave  equation,  is  statistically  homo¬ 
geneous,  at  least,  in  the  sense  of  second-order  co¬ 
herence  theory. 

It  is  also  apparent  from  eq s.  (2.8)  and  (2.14)  that, 
when  the  coherence  function  of  the  wavefteld  I'(/(,r) 
is  a  function  of  the  difference  variable  r  only,  the 
specific  intensity  I(R.s)  is  independent  of  position, 
i.e., 

/(R.s)  =  f(s).  (2.19) 

In  that  case  the  free  space  radiative  energy  transfer 
equation 

s-~Kl(R,s)  =  0  (2.20) 

is  trivially  satisfied. 

Using  these  results,  it  may  be  shown  that  the  rep¬ 
resentations  eq.  (2.17)  or  eq.  (2.18)  cannot  be  adopt¬ 
ed  to  derive  an  equation  of  radiative  transfer  with 
non-vanishing  extinction  and  scattering  coefficients 
in  statistically  homogeneous  media. 

In  this  paper  we  consider  wave  propagation  in  a 
statistically  inhomogeneous  medium  and  show  that 
the  representation  of  eq.  (2.17)  for  the  spectral  den¬ 
sity  f(R,K)  may  be  used  when  certain  assumptions 
are  made  concerning  the  incident  field  and  the  ran¬ 
domly  fluctuating  medium. 


3.  Wave  propagation  in  a  weakly  randomly  fluctuating 
medium 

The  propagation  of  the  space -dependent  part  :i  r\ 
of  the  scalar  wave  in  a  medium  with  random  fluctua¬ 
tions  of  the  dielectric  constant  is  governed  in  the 
stochastic  scalar  wave  equation 

+  A^ctnbVI  =  0.  (3.1) 

where  k0  =  u>:e  is  the  wave  number  of  the  field  in 
free  space  and  c(r)  is  the  dielectric  constant.  If  the 
fluctuating  dielectric  constant  exhibits  only  small 
departures  from  its  mean  value,  which  for  simplicity, 
is  assumed  to  be  unity,  then  a  representation  in  terms 
of  a  parameter  q  may  be  adopted  ; : 

eV)  =  1  +  qe,  (r).  (3.2) 

where 

q<(A0/,)-:.  (3.3) 

with  /]  denoting  the  smallest  correlation  distance  of 
the  dielectric  constant  fluctuation  e.  In.  Since  we  as¬ 
sumed  (e(r)>  =  1 .  we  have  (e,  (r)>  =  0.  Moreover  die 
random  function  e |  (z)  is  assumed  to  be  a  gaussian 
random  variable. 

On  substituting  from  eq.  (3.2)  in  eq.  (3.1 )  we  ob¬ 
tain 

(72  +  A.’oMr)  =  -k^pe^niir).  (3.4) 

If  we  recall  the  definition  (eq.  (2.1 ))  of  the  second- 
order  coherence  function  of  the  wavefteld  IV  . r - 1, 
it  follows  from  eq.  (3.4)  that  IV)  .z->  I  satisfies  the 
equation 

(V2  -  VjHV,  ,r2)  =  -*2q[<e,(/’1  )iV,  )r*(r, )' 

-  (e^/yliV,  )i'*(r,T]  ■  (3.5) 

Let  us  assume  that  the  field  variable  rtr)  can  be 
expressed  in  a  perturbation  expansion  [  1  '.IS] 

fV)  =  t'flfr)  +  qr>i  (rl  +  ()(q2  ).  (3.6) 

Upon  inserting  from  eq.  (3.6)  into  eq.  (3.4)  and 
equating  to  zero  die  coefficients  of  the  zeroth  and 
of  the  first  powers  of  q.  we  obtain 

(72  +k*)v0{r)  =  Q%  (57) 


Wc  assume  max.  <  t "  1 .  In  that  case  u  reprove  .he 
strength  of  the  dielectric  constant  fluctuations. 
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(7“  +  A'jj)p,(r)  =  -A3e,(r)t>0(/').  (3.8) 

The  solution  of  eq.  (3.8)  can  be  expressed  as 


n,(r)  =  -kljG0(r-p)ei(p)v0(p)dip.  (3.9) 


where 


G0(r-p)=- 


J 

4rr 


exp  (iA:0|/-  -p I) 
I r  -pi 


(3.10) 


is  the  (outgoing)  free  space  Green’s  function  of  the 
Helmholtz  equation  with  wave  number  A(1. 

We  now  expand  the  terms  on  the  right  hand  side 
of  eq.  (3.5)  to  the  lowest  non-vanishing  order  in  p. 
Since  <e,(r)>  =  0.  it  follows  from  that  perturbation 
expansion  (eq.  (3.6)),  (neglecting  the  terms  of  order 
p3),  that 


cj  -  722)r(r1.r2)  =  Aj/[G0(r1  -p)f0(p.r2) 

+  G*(r2  -p)r0(Vp)j  5(r,,p)d3p 
-k\  /[Gq(a,  -  p)f0(p.r2) 

+  Gj(r,  -  p)  r0(r, ,/»)]  5(p,r2)d3p,  (3.11) 


where 

B(r1,r2)  =  B(3(r1  +  /-,),/■  (  -  r2)  =  p3^,^,^,^)). 

(3.12) 

ro(r,>r2)=ro(|(r1+r2)(r,  -r2) 

-WlKW-  (3J3) 

We  now  change  to  the  variables  R  and  r  which  are 
defined  by  eq.  (2.2).  In  terms  of  the  new  variables, 
eq.  (3.11)  can  be  shown  after  some  rearrangements, 
to  take  the  form 

2VR-Vrr(«,r) 

=  k*  J[G0(p)B(R  +  \r-  ip.PH’o**  "  5 PS-P) 

-  G*(p)B(R  -  \r  - \p,p)\'0(R  -  \p,r  +p)]  d3p 

-  k*  J [G0(p)B(R  -5P.r-p)f0(«  - \p.r  -  p) 

-G*(p)B(R  —  {p,r  +  p)  rQ(f?  -  jp.r+p)]  d3p. 

(3.14) 


4.  Derivation  of  the  equation  of  radiative  energy 
transfer 

We  will  first  make  certain  assumptions  concerning 
the  incident  wavefield  v{)(r)  and  the  dielectric  con¬ 
stant  fluctuation  C|(r)  which  simplify  eq.  (3.14)  con¬ 
siderably.  We  then  proceed  to  derive  the  equation  of 
radiative  transfer  satisfied  by  the  specific  intensity. 

Let  us  assume  that  the  incident  field  is  statistical¬ 
ly  homogeneous,  at  least  in  the  sense  of  the  second- 
order  correlation  theory.  The  coherence  function 
r(,(/J.r)  then  depends  only  on  the  difference  variable 
r.  i.e., 

T()(7f  ,r)  =  r, ,(/■)•  (4.D 

Next  we  assume  the  random  medium  to  be  quasi- 
homogeneous,  i.e..  the  coherence  function  B(R,r ) 
of  the  dielectric  constant  fluctuation  is  assumed  to 
be  a  "slow"  function  of  R  and  a  "fast  "  function  of 
r.  If  /.  denotes  the  distance  over  which  B(R.r)  re¬ 
mains  sensibly  constant  with  respect  to  its  first  ar¬ 
gument.  then,  in  view  of  the  quasihomogeneity  of 
the  medium. 

L>lv  (4.2a) 

where  /,  is,  as  before,  the  smallest  correlation  dis¬ 
tance  of  the  fluctuations  in  the  dielectric  constant. 
Furthermore  we  assume  L  to  be  much  larger  than  the 
correlation  distance  /2  and  the  wavelength  2tt /A0  of 
the  incident  field,  i.e.. 

/,>/,.  r^2r:k0.  (4.2b.c) 

The*  av-v-mptions  are  satisfied  in  many  problems 
of  practical  interest.  A  frequent  model  to  describe 
a  random  medium  e.g.,  the  earth's  atmosphere  as¬ 
sumes  quasihomogeneity.  In  many  cases  of  optical 
propagation  through  atmosphere  the  conditions  (eq. 
(4.2a))  and  (eq.  (4.2c))  are  easily  satisfied.  The  as¬ 
sumptions  (eq.  (4.1))  and  (eq.  (4.2b))  are  obeyed  by 
fields  of  thermal  origin  e.g..  the  blackbody  radiation 
which  is  both  statistically  homogeneous  and  isotropic 
and  whose  correlation  extends  over  distances  of  the 
order  of  the  wavelength  associated  with  the  frequen¬ 
cy  of  the  radiation  at  which  the  energy  transfer  is 
being  considered. 

In  view  of  the  assumptions  (eq.  (4.1))  and  (eq. 
(4.2)),  we  can  simplify  the  various  terms  on  the  right 
hand  side  of  eq.  (3.14).  Since  the  incident  field  is 
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considered  statistically  homogeneous,  the  first  term 
can  be  written  as 

j  (;n(p)mR  +  \r  \p,r p)d3/> 

=  f G^p)B{R  +  kp.p)V0(r  -p)d3p.  (4.3) 

In  view  of  our  assumption  (eq.  (4.2b)),  B(R.r)  is  con¬ 
sidered  a  “slow"  function  of  R  and  1 (/■ )  is  consider¬ 
ed  a  "fast"  function  of  r.  Under  these  circumstances, 
we  can  replace  B(R  +  \r  -  Ip.p)  by  B{R,p)  on  the 
right  hand  side  of  eq.  (4.3)  without  introducing  an 
appreciable  error.  We  then  obtain 

[g0(p)B(R  +  Jr  -  J p,p)l'n(R  -  \p.r  p)dJp 

-J  G0(p)B(R,p)\'{](r  -  p)d3p.  (4.4) 

Similarly, 

J G*(p)B(R  -  i|-  -  ^p,p)l'0(K  -  \p,r  +  p)d3p 

-  j G*{p)B(R,p)\'0(r  +  p)d3p.  (4.5) 

Tlie  third  and  fourth  terms  on  the  right  hand  side  of 
eq.  (3.14)  can  also  be  simplified  by  using  eqs.  (4.1) 
and  (4.2b),  and  one  obtains 

J G0(p)B(R  -  {p.r -p)['0(R  -  \p.r  -p)d3p 

-  jcu(p)B{R,r  -  p)\'0yr  -  p)Aip.  (4.6) 

fc*(p)B(R  -  \p.r  +p)rn(R  -ip,r+p)d3p 

*  jG*(p)B(R.r+p)r0(r+p)d3p.  (4.7) 

On  substituting  from  eqs.  (4.4)-(4.7)  in  eq.  (3.14) 
we  obtain 

•  7,r(*,r)  =  k*  JlG0(p )  l'0(r -p) 

-G*(p)  r0(r+p)]  B(R,p)d?p 
-k*J  [G0(p)B(R.r -p)\'0(r-p) 

-  G*(p)B(R,r  +  p)  T0(r  +  p)|  d3p.  (4.8) 

If,  on  the  right-hand  side,  r0(r  +p)  is  replaced  by 
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\'{R.r  z  pi  then,  in  view  of  eq.  (3.01.  an  error  of  the 
order  p4  is  introduced  which  is  neglected  in  the  pres¬ 
ent  lowest  non-vanishing  order  perturbative  analysis 
[17],  Thus,  to  the  second  order  inp. 

27^-  7rr(*.r)  =  *2/[G„(pir(*.f  pi 

G*(p) \'{R.r  +p)]  B(R. p\d~'p 

-  k*f  [G0(p)B(R.r  p)l  (R.r-  p ) 

-G*(p)/i(/?.r+p)l(/t.r+p)]  d3p.  (4.9) 

Next  we  take  the  Fourier  inverse  of  eq.  (2.S)  to 
express  l  (R.r)  in  terms  of  the  spectral  density/)/?. 

AT),  i.e.. 

F(/?.r)=  //(/?, A)exp(iAV)d3A.  (4.10) 

Moreover  we  define  the  spectral  density  '1>(/?.  A)  of 
the  function  B(R,r)  by  means  of  Fourier  decomposi¬ 
tion  with  respect  to  r. 

B(R,r )  -  J  'I>(/?,A)exp(iAT)d3A.  (4.1 1 ) 

It  can  be  shown  that  <!>(/?.  A)  is  a  real  function.  In¬ 
serting  eq.  (4.10)  in  eq.  (4.9)  and  using  eq.  (4.1 1 ). 
we  obtain  after  some  rearrangement 

K-1Rf(R.K)  =  -JjrA*  jj'hlA.A  -  A-0s’)df2s.|/(/?A') 

+  j  jtAq  S(A  -  A0)  f -!-(/?. A  -  K')r\R.K')d:'k". 

(4.12) 

This  integro-differential  equation  for  HR.  A)  can 
be  recast  into  the  following  integral  equation  1 : : 

f(R.K)=  .(j((A)  exp  J  -  J  dx"  A(R  -  x"s.  A)1 

'  o  ■ 

+  6(A-k0)  J  d x'B(R  .v  s, A) 

0 

X  exp  J  f  dx",l(R  -.v'i.A)'.  (4.13) 

'  o  I 

* 3  The  derivation  of  eq.  (4.16)  can  be  carried  out  on  the 
same  lines  as  in  Ch.  2  (Sec.  6)  of  ref.  [  2 ) . 


OPTICS  COMMUNICATIONS 


319 


OPTICS  C'OMMI  NK  A 1  IONS 


1  June  1 MH 1 


Volume  3?.  number  5 

where 

/..( K)-  — -jl(r)exp(  \K-r)il}r.  (4.14) 

CrrE 

nkl)  C 

A{R.K)  =  j'HR.K  -k{)s  (4.1M 

~ko  - 

«(K.A)  =  —  I  MR. A  K')/(R.R‘)dK'.  (4.16) 

—A  • 

Since  the  incident  field  is  assumed  to  be  statistically 
homogeneous,  it  can  be  shown  that  (ret'.  [12| .  appen¬ 
dix  Bt 

f0(X~)  =  /|t(s)iS(A'  -  kf))/k~t.  (4.17) 

On  substituting  eq.  (4.17)  in  eq.  (4.13).  it  becomes 
obvious  that  J(R.K)  is  of  the  form 

f\R.K)  =  l(R.s)6(K  k  „)/*;•.  (4.18) 

This  relationship  between  the  spectral  density  f\R. 

K)  and  the  specific  intensity  HR.s).  which  is  a  na¬ 
tural  consequence  of  the  assumptions  made  in  this 
and  the  previous  sections,  has  been  shown,  in  sec¬ 
tion  2.  to  satisfy  the  "space  density'1  relation  (eq. 

(2.1 1))  and  the  energy  flux  relation  (eq.  (2.12))  of 
the  phenomenological  theory. 

If  this  expression  for  f(R.K)  is  now  substituted 
in  eq.  (4.12),  then  after  integration  over  K' .  [d3A’ 

=  K'-dK' | .  we  obtain  the  following  equation: 

s-7rI(R.s)=  -a(R.s)/(R.s) 

+  f  0(R.s.s')HR,s')  dSl4,,  (4.19) 

where 

a(R.s)='ink*J>liR,'k0{s  -s')) (4.20) 

HiR.s.s')  =  \rtk40MR.k{)(s  -  s’)).  (4.21) 

Eq.  (4.19)  is  of  the  same  form  as  the  equation  of 
radiative  energy  transfer  of  the  phenomenological 
theory  (2| .  Our  formulae  (eq.  (4.20))  and  (eq.  (4.21 )) 
relate  the  extinction  coefficient  and  the  scattering 
coefficient  0  to  the  stochastic  '-hiral-fro<-t;LS  of  the 
medium. 

It  is  apparent  from  eqs.  (4.20)  and  (4.21 )  that  a 


and  3  are  connected  by  the  relation 

c4R.s)=f 0(R.s.s')d<ls.  (4.22) 

This  equation,  which  is  a  consequence  of  the  conser¬ 
vation  of  energy,  implies  that  the  extinction  of  a 
wave  propagating  along  the  direction  s  is  due  to  the 
scattering  in  all  other  directions. 

It  is  a  pleasure  to  acknowledge  the  encouragement 
and  assistance  of  Professor  Emil  Wolf  during  the 
course  of  this  work.  The  author  is  also  grateful  to 
Professor  M.O.  Scully  for  hrs  support. 
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The  steady  dale  quantum  statistical  properties  nl  a  two-photon  laser  are  presented  tor  the  case  when  the  emus  loss 
mechanism  is  simulated  hv  a  simile-photon  process. 


I.  Introduction.  The  success  of  the  quantum  theory  of  lasers  involving  a  single  photon  emission  per  atomic 
transition  [  I  ]  has  generated  a  great  deal  of  interest  in  the  possibility  of  achieving  laser  action  involving  the  stimu¬ 
lated  emission  of  two  or  more  photons  in  a  single  atomic  decay.  The  two-photon  laser  has  the  prospect  of  achiev¬ 
ing  high  beam  intensities  since  the  strength  of  the  coupling  between  the  laser-active  atoms  and  the  light  field  is 
proportional  to  the  light  intensity  in  contrast  to  the  square  root  of  the  light  intensity  in  the  one-photon  laser. 

Several  authors  [2-7 1  have  studied  the  photon  statistics  of  a  two-photon  laser  on  the  basis  of  a  model  similar 
to  the  laser  model  of  Scully  and  Lamb  [  1 1.  The  model  consists  of  a  coupled  system  of  a  field  and  identical  two- 
level  atoms.  The  lasing  levels  of  the  atoms  are  assumed  to  have  the  same  parity  under  the  usual  dipole  approxima¬ 
tion  for  the  two-photon  transition  to  take  place.  The  atom— field  interaction  for  the  two-photon  emission  process 
may  be  described  by  the  effective  hantiltonian  [8]: 

H\  =g2<o+J2  +  o  ~a*2) ,  (1) 

where  g,  is  the  coupling  constant  given  by 

,?2  =  ju  1 2|/:1  -  -  . 

ju  ]  is  the  matrix  element  for  the  two-photon  transition,  and  |/:'|  2  has  the  dimension  of  field  intensity,  the  atom 
operators  o  and  of  are  defined  by 


,0  1  s 

,0  0. 

• 

H 

'()  o' 

'  1  0' 

and  a  .a*  are  the  annihilation  and  creation  operators  of  the  field.  The  cavity  losses  due  to  the  transmission  of 
laser  light  through  the  end  mirror  are  simulated  by  another  set  of  two-level  atoms  which  are  pumped  in  the  lower 
level  at  a  constant  rate  and  which  can  make  a  transition  to  the  upper  level  by  absorbing  laser  photons. 

In  most  of  the  studies  on  the  quantum  statistical  properties  of  the  optical  field  ol  a  two-photon  laser  1 2.5,t>] . 
it  is  assumed  that  the  end  mirror  of  the  cavity  is  transmitting  at  the  sum  frequence  2 co.  i.e..  the  two-level  atoms 
which  simulate  the  cavity  losses  absorb  two  photons  to  make  a  transition  from  the  lower  level  to  the  upper  level. 
Tlie  restriction  to  a  two-photon  loss  mechanism  is  made  in  order  to  retain  the  property  of  detailed  balance.  This 
model  for  the  cavity  losses  is  unrealistic  since  the  photons  do  not  escape  in  pairs.  On  the  basis  ol  this  model,  it 
is  predicted  that,  in  single-  and  two-mode  two-photon  lasers,  the  photon  distribution  'mictions  are  narrower  than 
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in  the  single-photon  laser. 

It  ts  of  considerable  interest  to  investigate  the  quantum  statistical  properties  of  a  two-photon  laser  when  the 
cavity  loss  mechanism  is  simulated  by  a  single-photon  process.  Some  authors  have  approached  this  problem  ustne 
perturbation  methods  [3,7|.  In  this  paper,  we  present  an  exact  solution  ol  the  master  equation  that  describes  a 
single-mode  two-photon  laser  with  a  single-photon  loss  mechanism.  We  als<  derive  expressions  tor  the  mean  pho¬ 
ton  number  and  the  relative  fluctuations  of  the  number  of  photons  hii-1'  above  threshold. 


2.  Photon  distribution  Junction.  We  start  with  the  equation  of  motion  for  the  photon  distnbution  function 
p(n)  in  a  single-mode  two-photon  laser.  The  master  equation  can  be  obtained  in  a  straight  forward  manner  h>  fol¬ 
lowing  the  same  method  as  that  of  the  Scully-Lamb  theory  of  a  single-photon  laser  1 1  ].  The  resulting  equation  is 


<}/>(«> 

dt 


.1-lu  +  1  Hn  f  2>  A  iinn  1  » 

I +(/!,  .1  •>)(// +  1  )(u  +  2)  ^  1  +(/#_>  -lifutn  \)  ^" 


2)  +  C(n  +  1  )p{n  +  1 )  Oipin  I . 


(2) 


where,  similar  to  the  Scully  -  Lamb  theory,  .1 1  and  B-,  are  the  gain  and  saturation  patameters  tor  the  two-photon 
gain  mechanism  and  C  is  the  loss  parameter  which  is  related  to  the  laser  frequency  ui  and  the  (J  ot  the  cavity  by 
the  following  equation: 


C  =  ijo/Q  . 


(2 1 


In  eq,  (2),  the  term  proportional  to  p(n  -  2)  represents  the  two-photon  emission  process,  whereas  the  term  pro¬ 
portional  to  p(n  +  1 )  represents  the  single-photon  absorption  process.  It  is  therefore  evident  that  the  detailed 
balance  condition  is  not  obeyed  in  the  present  case.  This  makes  it  more  complicated  to  solve  eq.  (2) even  in  the 
steady  state. 

In  an  earlier  paper  (9)  on  the  effect  of  cooperative  atonuc  interactions  on  photon  statistics  in  a  single-mode 
laser,  we  encountered  a  master  equation  which  also  included  two-photon  emission  and  single-photon  absorption 
processes.  We  solved  this  equation  in  steady  state  using  a  matrix  approach.  Hie  same  approach  can  be  applied  to 
the  present  problem. 

In  the  steady  state  [dp(n)/dt  =  0],  we  obtain  from  eq.  (2). 
anp(n)  +  bn_2P(n  -  2)  +  c„+ ,p(/t  +  1)  =  0  ,  (41 

where 


__  (n  +  i)(n  +  2)_  =  <m +!)<«  + 2) 

a"  f  +  (*2//t ,)(«  +  ' n(/t  +  2")  •  n  1  +<S,/,l:)(n  +  l)(;t  +  2)  ' 

It  is  clear  from  eqs.  (5a-c)  that 
an  +  bn  +  cn  =  0  • 

In  matrix  notation,  eq.  (4)  can  be  rewritten  as 


('a.b.cl 


((’» 


ja°  c'  o\  }H0) 

0  1 /,  c-2  ,  i  i»( I ) 

>b„  o\  ••  j  ; 


It  can  be  shown,  by  the  method  of  induction,  that  the  solution  of  eq.  ( ?)  is  given  by  the  following  equation. 


226 


\  v'lumc  ho.\.  number  4 


PinSK  S  I  I  1  11  HS 


s  hcicmtvr  l‘*S" 


(  DVUM 

/><«»  =  Mon. 

<1  A  ...  <  t* 


where 


\Hii )  =  ile  I 


“n  <  i 

0  u. 


1  n  I 


0  a 


>i  1 


<S> 


<lM 


In  eq.(S) p(0)  is  detei  mined  irom  the  normalization  condition  )  =  I .  By  using  the  properties  of  the  de¬ 

terminant  and  i  elation  1 01.  the  determinant  Min )  can  he  shown  to  he  expressible  in  the  I  olio  wine  form: 


Min  I  =  del 


On  substituting  trom  eq.  ( 101  in  eq.  (is).  we  obtain  />(u )  as  the  loiiowrng  product  of  continued  tractions: 


(10) 


n 

pin)  =  p<()>  n  ;v 

r  I  •> 


f'r  :<■  r  1 


hr  i  + 


f’r  Mr  2 


Tins  expression  tor  pin )  together  with  the  expression  ol  /’„  and  <  „  |cl.  eqs  (5h,c))  completely 
photon  statistics  ot  a  two-plioton  laser 


(II) 


determines  the 


s’.  I’h, ‘inn  number  Ilia  riuihnm  huh  tlm  shnU.  We  now  consider  the  photon  distnhution  lutistion  ot  the 

two-photon  laser,  when  the  laset  operates  Inch  above  threshold,  a  regime  in  which  the  photon  distribution  tunc- 
tton  ot  a  single-photon  laser  is  described  In  a  I’oisson  distribution,  l  or  a  I’oisson  distnhution.  the  tluctuation  in 
the  number  ot  photons  is  unity  .  i.e.. 

■r'-v-- 

boi  a  two-photon  laser,  opeiat me  high  above  threshold,  we  have  n  ■  (  I  s  Hi )'  -.  Tinier  this  condition  eq.  ( 
can  be  simplified  considerably  by  neglecting  I  in  the  denominators  ol  the  lust  and  second  terms  ol  the  left  hand 
side.  We  then  obtain 

dp(«  I  dr  -  (.  I  %  Ws  t/nn)  +  t.-ly  H -)/>(».  >  '  i  in  t  I  )pin  I )  Or/>(»r).  (Id) 

Tins  equation  can  easily  be  solved,  in  steady  state,  using  the  following  generation  function: 
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£(.t)=  Tj  x"p(n) . 

(14) 

H  =  () 

The  generating  function  g(.v )  obeys  the  following  equation: 

-  (.-15/B2C)(1  +  .v\?  =  0. 

(15l 

This  equation  can  be  integrated  and  the  resulting  equation  is 
g(.v)  =  expKA^flsOl-V'  +  lv  -  3)|  . 

It  is  evident  from  eq.  ( 14)  that 

(16) 

<„>  =  ?'  ,  <„2>«/-d2*  4^)  . 

d-v  \dr2  dx/A.=  1 

(17a,  b) 

It  follows,  on  substituting  from  eq.  (16)  in  eqs.  (17a,b)  that 
(n)=2A\/B:C  ,  (An)2/(«>  =  ^  . 

(IHa.bl 

According  to  eq,  (18b),  the  photon  distribution  function  for  a  two-photon  laser  is  wider  than  the  Poisson  distri¬ 
bution.  This  result  is  in  agreement  with  that  of  Golubev  [7f  that  the  relative  fluctuations  of  the  number  of  pho¬ 
tons  in  a  two-photon  laser  cannot  be  smaller  than  in  coherent  emission. 


It  is  a  pleasure  to  thank  Professor  M.O.  Scully  for  his  support. 
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A  density  matrix  method  is  used  to  obtain  an  exact  solution  for  the  reduced  density  matrix  of  the 
field  in  an  arbitrary  multimode  m-photon  absorption  process.  The  results  of  some  earlier 
specialized  studies  of  photon  statistics  in  multiphoton  absorption  process  can  be  recovered  from 
this  solution. 


It  is  well  known  that  nonlinear  interaction  of  light  with 
matter  changes  the  quantum  statistical  properties  of  the 
light  field.1  The  changes  of  the  photon  statistics  depend  on 
the  interaction  process  and  on  the  initial  conditions  of  the 
field. 

In  this  paper  we  consider  a  multimode  m-photon  ab¬ 
sorption  process.  The  multiphoton  absorption  processes 
have  recieved  a  great  deal  of  attention  in  recent  years  due  to 
the  possibility  of  producing  a  radiation  field  which  shows 
photon  antibunching.  The  photon  statistics  of  two-photon 
absorption  process  was  studied  using  a  generating  function 
approach  and  exact  expressions  for  the  photon  distribution 
function  were  given  in  single-mode"  '  and  two-mode"  pro¬ 
cesses.  The  off-diagonal  elements  of  the  density  matrix  of  the 
field  in  the  single-  and  two-mode  two-photon  absorption 
processes  were  given  by  Simaan  and  Loudon.7 "  Paul,  Mohr, 
and  Brunner'  studied  the  photon  statistics  of  m-photon  ab¬ 
sorption  process  on  the  basis  of  an  approximate  procedure. 
Recently,  an  exact  analytic  solution  of  the  master  equation 
describing  single-mode  m -photon  absorption  has  been  ob¬ 
tained  by  Zubairy  and  Yeh1"  using  a  density  matrix  ap¬ 
proach  and  by  Voigt,  Bandilla.  and  Ritze"  using  a  Laplace 
transform  method. 

In  the  present  paper  we  extend  these  results  to  arbitrary 
multimode  absorption  process.  We  begin  by  considering  the 
master  equation  which  describes  a  multimode  m-photon 
process.  We  present  the  exact  solution  of  the  reduced  density 
matrix  of  the  field  using  the  matrix  approach  Due  to  the 
general  nature  of  our  problem  the  results  of  the  earlier  stud¬ 
ies  can  be  recovered  in  the  appropriate  limits. 

We  consider  a  coupled  system  of  a  field  and  A  noninter¬ 
acting  two-level  atoms  in  their  ground  state.  The  number  of 


atoms  ,V  in  the  lower  level  are  assumed  to  be  maintained 
constant  by  some  external  influence.  We  assume  that  the 
atoms  make  a  transition  from  the  lower  level  to  the  upper 
level  by  absorbing  m  photons;  jc,  photons  in  mode  i  (i  =  1.2. 
....  /;  Km).  We  then  obtain 


The  single-mode  situation  corresponds  to  the  case  when  one 
of  the  x,’s  such  as  x,  is  equal  to  m  and  all  the  rest  are  zero. 
The  reduced  density  matrixp,  of  the  field  can  then  be  shown 
to  satisfy  the  following  equation  of  motion  using  the  stan¬ 
dard  perturbation  techniques 

x  (,n,a:  • ) + a  (,n,«;  ■ )( n,---  )]■  ,J» 

where/?  ’  is  the  absorption  coefficient  for  w-photon  absorp¬ 
tion,  and  d'.  a,  are  the  photon  creation  and  destruction  oper¬ 
ators  of  the  /th  mode,  respectively.  In  Eq.  |2|  the  saturation  of 
the  absorbing  atoms  is  neglected. 

We  denote  the  state  in  which  there  are  n  photons  in  the 
/th  mode  |»  =  1,  2, ....  / )  by  zr,,  n  ..n  ).  The  equation  of 
motion  |2|  for/5,,  can  be  translated  into  an  equation  for  the 
matrix  element 

pin, . n,\  n ,  +  A . . n,  +  A,.  r| 

=  (n, . n,\pt  ,n,  +  A, . n,  +  A,),  |3) 

of  the  reduced  density  matrix  by  evaluating  the  matrix  ele¬ 
ment  of  each  term  in  Eq.  |2)  between  the  appropriate  Fock 
states.  The  resulting  equation  is 


dp(n . .  n,  4-  A , . n ,  -I-  A,;r) 

dr 


where  r  -2P{m)<  and 


=  a(n . .  ft,;  n,  +  A . n ,  -t-  K,)p(n ,, ....  ft,;  n,  -t-  A,, ....  n,  +  A',;  r) 

+  />(«,  +  x . . «,  +  x,;n,  -f  x,  +■  A, . n,  -t-  x,  +  A',) 

Xp(n |  +  x, . n,  4-  x,;  n,  4-  x,  +  A, . n,  +  x,  4-  A,;  r), 


. . . . 


(4) 

(5a) 
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Following  the  method  employed  in  Ref.  10,  we  first  express  Eq.  (4)  in  the  matrix  notation 


~~p  ,( r|  =  A#  tp  ;(r|, 

dr 


where 


P.. 


p(ti.  i,;  i  i  4-  A^i . *;  4-  Kp,  r| 

p\i ,  4-  x , . i,  4-  x,;  t 1  4-  AT ,  -t-  x,. l;  4-  /C,  4-  x, ;  x) 

y>|i,  4-  nx . i.  4-  «x,;  i,  4-  K,  4-  nx„  ...  /,  +  AT,  +  nx,,;  r) 


m,K) 


M,  ,  = 


r 


l5bl 


|6I 


1 8) 


In  Eq.  (8),  the  elements  m,k  are  given  by  the  following 
expressions 

m„  =  <*(/',  4-yx, . i,  4-yx,; 

i,  +  K,  -t-Jx, . i,  +  K,  4-yx,),  |9a| 

mn  .1  =  b  If  I  j  4-  Dx, . i,  +  I  j  +  Ux,; 

i\  4-  K ,  -t- 1  j  ■¥  1  Lx, . i /  4 -  K,  4-  (y  4-  1  Lx,- 1,  l%l 

fory  =0.  1.2.-. 

The  solution  of  Eq.  (6|  or  Eq  |4|  can  be  expressed  in  the 

form 

p[i,  4-  nx,. ....  i,  4-  nx,;  /,  4-  A,  4-  nx„  ....  i,  4-  K,  4-  nx,;r| 

=  £  X  . '/  +1*ii 

>  l>  V  <* 

Xi,  +  K,+qx, . i,  +  K,  4-  ?x,;  0),  (10( 

where  A,  ,,,  ,  are  the  eigenvalues  of  the  matrix  A/,..  , .  i.e.. 
they  satisfy  the  equation, 

det[Af,  ,  -  X/  ]  =  0.  (Ill 

(/being  the  unit  matnxl,  a’'”  is  the  nth  element  of  the  right 
eigenstate  of  M,  ,  corresponding  to  the  eigenvalue /i, ,  ,  . 


and /S’’"  is  the^th  element  of  the  left  eigenstate  of  M 
corresponding  to  the  eigenvalue / ,  .  The  matrix  ele- 
menlsa,/’  ''and/?]/1  can  beshown  uiohoy  the  following 
recursion  relations 

m„  i*  i  <*!/"  i'’’ 4- m„  il2ai 

m,  ,v0'v'" +  m„0l"  '  1 ' -hi 

By  solving  Eq.  (11)  and  by  iterating  the  recursion  relations 
( 1 2a)  and  (12b).  we  obtain 


(l.M 

!  FT  (  ‘ 

1  n-  t 

=  1 

|-  V.  i\m„  -  m. 

1,  N  V  Ji 

/  U4) 

L.0;  n  ->  s. 

ft  (  m-  " 

)  on 

K 

—  \ 

,  M.  i\m„  -  m„ 

If  we  let  n 

=  0,/ 

.0;  q  <  s. 

,  =  n  ,  and  A.'  -  n  ,' 

-  n  (y  -  |,  2  7  )in  Ei 

(10),  then,  on  substituting  from  Eqs.  (9a).  (9b).  and  ( 13)- 
(15).  it  follows 


pin, . r) 


-l  1 


II?  ,b(n,  4- ex, . n  4- rx,;n;  4 -  rx, . n‘,  -t  rx,) 


,)[a(n,  4-  rx, . n,+sx,;n I  4-  sx . »;  4- xx,)  -  afn,  4- «t„  ....  n,  4- rx,;  n;  4- rx,. ....  «;  4- rx,)  j 


Xe 


*■  »«,.  *ij  ♦  *4,)r 


pin,  +qx, . n,  4 -qx,\n\  4 -qx,,...n;  4-<?x,.0). 


(lb) 


This  equation,  combined  with  the  expressions  fora’s  and  b ’s 
[cf.  Eqs.  (5a)  and  (5b)],  completely  determines  the  time  evo¬ 
lution  of  the  density  matrix.  The  photon  distribution  func¬ 
tion  pin . r)  =  (n . n,\pf  ....  n,)  can  be  deter¬ 

mined  from  Eq.  (16)  by  putting  n =  n,  (y  =  1, ..../). 


I 

We  now  show  how  some  of  the  earlier  specialized  re¬ 
sults  regarding  multiphoton  absorption  processes  can  be  ob¬ 
tained  from  Eq.  (lb). 

In  the  case  of  single-mode  m-photon  absorption  pro¬ 
cess,  one  of  thex.’s  (say  x,)  is  equal  tom  and  all  the  rest  arc 
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--l/_»a! — -  - ^ - \  (.7a) 

2  \(n i  —  m )*  («',  -  m)\J 


zero.  We  then  obtain,  from  Eqs.  (5a)  and  (5b).  that 
«(«,;  n\ )  = 
pin,',  n\\  r) 

=  S  I 


/  *,!«;!  V  - 

*(»,;»;)  =  ( - — - 

V»,  -«)!(«;  -  w)!/ 


[17b) 


If  follows  from  Eq.  (16)  that 


Ii;  tb(n,  +  mr,  n\  +  mr) 


a(n,  t  ms.n't  »  "i 


,11;  ,,[<j(/i|  +  ms;  n\  +  ms)  —  a(nx  +  mr\  n\  +  mr)] 


’  """pin |  +  mq\  n[  +  mq\  0). 


(18) 


After  making  some  rearrangements,  this  equation  can  be  shown  to  be  identical  to  Eq.  (31)  in  Ref.  10  (see  also  Ref.  1 1  for  the 
case  «,  =  n\ ).  We  have  discussed  some  aspects  of  photon  statistics  in  single-mode  m-photon  abosorption  process,  such  as 
photon  antibunching,  in  that  paper. 

As  another  example,  we  consider  m-mode  m-photon  absorption  process.  In  this  case,  /  =  m  and  x  =  l(i  =  1,2 . m). 

The  expressions  for  a  and  b  in  Eqs.  (5a)  and  (5b)  then  become 

ain . .  n', . «,'„)=  -([«,. ...  nm  n'm  ].  (19a) 

bin . .  n\ . n'm)  =  (n,  (19b) 

Moreover,  from  Eq.  (16),  we  obtain  the  following  solution  for  the  density  matrix 


pin . .  n\ . r) 

=  s  i 


,  b  (n  i  +  r . nm  +  r,  n\  +  r . n„  +  r) 


H;  „[u(ti|  +  s, . nm  +s;  n\  +  s . n'„  +  s)  -  a(n,  4-  r, ....  nm  +  r;  n\  +  r, ....  n‘m  +  r] 


X-e*’-  '  ’•  ‘  '"pin,  +  $ . nm  +q;  n\  +  q . n’m  +  q\  0). 


Simaan  and  Loudon'' K  have  discussed  the  case  m  =  2  in 
some  detail.  Equation  (20)  together  with  Eqs.  (19a)  and 
(19b),  after  some  rearrangement,  can  be  shown  to  reduce  to 
their  results.  For  a  discussion  of  the  quantum  statistical 
properties  of  the  double-beam  two-photon  absorption  pro¬ 
cess,  we  refer  the  reader  to  Ref.  6. 

In  conclusion,  we  have  obtained  an  exact  solution  of  the 
master  equation  that  describes  an  arbitrary  multimode  mul- 
tiphoton  process.  We  have  shown  via  two  examples  how  the 
results  of  the  earlier  specialized  studies  can  be  recovered 
from  this  general  solution. 


(20) 
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Intensity  correlations  in  a  two-mode  laser  oscillating  on  two  coupled  transitions  have  been  studied  by  means  of  the 
FokkerPlanck  equation.  It  is  shown  that  they  can  be  expressed  in  terms  of  the  eigenvalues  and  eigenfunctions  of  a 
one-dimensional  Schrodinger-type  equation  It  is  found  that  near  threshold  intensity  correlations  cannot  be 
approximated  by  a  single  exponential.  For  moderately  large  excitations,  however,  a  single  exponential  dominates. 
Approximate  analytic  expressions  for  large  excitations  have  been  obtained  which  show  that  the  correlation  time 
increases  with  increasing  excitation. 


I.  INTRODUCTION 


Recently  a  quantum- mechanical  treatment  of 
laser  oscillations  on  two  coupled  atomic  transitions 
has  been  given.1  In  that  treatment  a  set  of  homo¬ 
geneously  broadened  three- level  atoms  interacting 
with  a  two- mode  electromagnetic  field  in  a  uni- 
directional-ring  configuration  was  considered. 

The  atomic- level  configuration  that  was  con¬ 
sidered  is  usually  referred  to  as  the  “A"  con¬ 
figuration  in  the  literature  on  three-level  atomic 
media.2  It  was  assumed  that  each  transition  sup¬ 
ported  a  distinct  mode  of  the  electromagnetic 
field.  A  master  equation  lor  the  photon  distribu¬ 


tion  function  was  derived  under  the  conditions 
of  perfect  resonance  between  the  atomic  transition 
and  the  corresponding  mode  frequencies.  This 
equation  was  solved  in  the  steady  state  and  the 
fluctuation  properties  of  the  optical  field  were 
discussed.  One  of  the  principal  conclusions  was 
that  the  intensity  fluctuations  do  not,  in  general, 
die  away  with  increasing  excitation.  Small  asym¬ 
metries  in  gain  may  change  the  statistical  proper¬ 
ties  of  the  two  modes  significantly.  These  effects 
are,  of  course,  the  result  of  mode  competition. 
The  question  that  presents  itself  now  is  how  does 
the  mode  competition  affect  the  correlations  and 
other  time-dependent  phenomena.  This  paper 
answers  this  question  partly. 


II.  TIME  DEPENDENT  FOKKER  PLANCK  EQUATION 

The  starting  point  of  our  discussion  is  the  master  equation  for  the  photon  distribution  function  /><",, //,) 
derived  in  Ref.  1,  viz., 

d  .  A(n,  + 1)  .  ,  A[n,+ 1)  .  . 

Jt /»(»!. »i>  +  +  A)(n,  *nz  +  2)P(n"  2 

+  C,  (n,  +  l)/>(n,  + 1 ,  n2)  +  C2(n2  + 1  )p(nx ,  n2  + 1 ) 

+ 1  +  (57 A)(n\  *n2  +  "j)  + 1  +  (B/A)(>,*+  "2 '  l' 

-Clnlp(nl,nt)~C2ntp{nl,n,) .  (D 


Here  />(«,, N1)  =  p(M1«ra;»,«a;f)  is  the  diagonal  element  of  the  density  matrix  operator  of  the  optical  field; 
n,,Hj  are  the  occupation  numbers  of  the  two  modes  and  A,  B,  and  C,  are  the  gain,  saturation,  and  the  loss 
coefficients,  respectively.1  This  equation  will  now  be  converted  into  an  equation  for  the  quasiprobability 
distribution  tunction  <l>  for  the  complex  field  amplitudes  by  using  the  coherent- state  representation  of  the 
density  matrix  of  the  electromagnetic  field.1  The  Fock-state  elements  of  the  density  matrix  operator  are 
related  to  4*  by 


f 


d3d,d%  , 


(2> 
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where  is  the  complex  field  amplitude  associated 
with  the  it h  mode  in  the  coherent  state  of  the 
field.  In  what  follows,  our  interest  will  be  con¬ 
fined  mainly  to  the  intensity  correlations  so  that 
we  will  not  be  concerned  with  the  dependence  of 
*l>  on  the  phases  p1  and  p2.  Such  dependence  would 
be  needed  for  a  discussion  of  the  amplitude  corre¬ 
lations  but  the  problem  in  that  case  turns  out  to  be 
rather  complex.  Now,  the  relevant  intensity  dis¬ 
tribution  function  P(/,,/2,f)  that  is  needed  to  cal¬ 
culate  intensity  correlations  can  be  projected  out 
of  4>  simply  by  noting  that  it  is  determined  solely 
by  the  diagonal  Fock-state  matrix  elements  of  the 
density  operator.  The  two  are  related  by 

P(« i«2;  «,«2; t)  t ) 

f  A  A 

(3) 


for  ic  or  p.  However,  if  we  take  recourse  to  the 
scaling  argument  a  closed  equation  for  I1  can  be 
obtained.  This  argument  works  as  follows:  since 
Bi  A  is  typically  -10"*  so  that  to  a  very  good  ap¬ 
proximation  the  derivative  terms  (B  'A)lt(i  s/,1 
can  be  ignored  compared  to  ( B ,  /!)/,  in  Eq.  (5b) 
which  then  is  easily  solved  to  give 

M=r  +(B/A)(ii*i2)  (6) 

On  substituting  Eq.  (6)  into  £q.  (5a)  we  obtain  the 
following  closed  equation  for  P(/,,f2,  M: 


dP_ 
dt  ' 


I  3  3  3  3 

_A\3A;‘  +  bIt  2~  dA/l  d/i 

*(c,a 


(7) 


We  also  introduce  an  auxiliary  function  by 
»('!,.  »,)  =  (l+^(«i+«a  +  2)) 

=  Jw(f„/2,<)e-,»-*^7-^-«Hlrf/2.  (4) 

«2. 

Using  Eq.  (4)  we  obtain  the  following  set  of  coupled 
differential  equations: 

(5a) 

and 

[1_  m(/i if^i ' 

(5b) 

These  equations  can  be  solved  exactly  in  the  steady 
state  and  the  corresponding  solution  was  discussed 
in  Ref.  1.  The  general  time- dependent  problem 
for  p(fltIt,l)  is  very  complicated  because  it  does 
not  seem  feasible  to  obtain  a  closed  equation  either 


Note  that  Eq.  (7)  has  the  formal  structure  of  a 
Fokker- Planck  equation  and  is  the  desired  equa¬ 
tion  that  will  form  the  basis  of  the  discussion  in 
the  rest  of  the  paper.  We  emphasize  that  the  ap¬ 
proximation  made  to  reach  at  Eq.  (7)  is  not  the 
weak-signal-limit  approximation  and  does  not 
restrict  the  validity  of  the  above  approach  to  the 
region  near  threshold. 

The  effects  of  mode  competition  are  more  inter¬ 
esting  when  both  the  modes  have  equal  losses.  We 
shall  therefore  consider  this  case  first  and  put 
C,  =C2=C  in  Eq-  (7).  The  general  problem  with 
unequal  losses  is  still  quite  complicated.  We  shall 
consider  this  problem  later.  With  equal  losses 
for  the  modes  we  now  look  for  the  general  solution 
of  the  form 

,  (8) 

Lm 

where  C  Lm  are  the  coefficients  to  be  determined 
by  the  boundary  conditions  and  \tm  is  the  eigen¬ 
value  of  the  differential  operator  on  the  right-hand 
side  of  Eq.  (7)  associated  with  the  eigenfunction 
/ Lm.  The  differential  operator  depends  on  two 
variables  /,  and  J2,  therefore  the  eigenvalues  and 
the  eigenfunctions  are  labeled  by  two  integers 
L,m.  The  eigenvalue  equation  is  obtained  from 
Eqs.  (7)  and  (8)  to  be 


f  Lm  ^ Lmf L 


In  order  to  solve  Eq.  (9)  we  first  try  a  substitution  of  the  form 


(9) 


(10) 

where  is  some  function  of  f,  and  I,  which  will  be  chosen  to  yield  a  self-adjoint  equation  for  g  Lm  and 

N  is  a  normalization  constant.  It  is  easy  to  show  that  the  choice 
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/■(/„/,>  =  -  m(i  +£(/,♦  (/,  +  /j  (/,  w/ 

yields  the  following  self-adjoint  equation  for  i*Lm: 


Notice  that  the  steady- state  solution  corresponds 
to  L  =  0  =  »i ,  A  tm  =  0  and  is  found  from  Eq.  (9)  after 
straightforward  integration  to  be 

/».Ui  ,/I)=Q-‘(l+J  (/,  +  /»)) 

xexp(i^(/l  +  /2)-g(;lw2r') .  (131 

If  we  choose  ,V  =  \  ij,  the  relation  (Eq.  (10))  can  be 
rewritten  for  later  convenience  in  the  form 

+  .•l)(ll+/l)J‘«Jft.(/lI/g).  (14) 

The  functions  may  be  chosen  to  be  an  ortho¬ 
normal  set.  In  addition  we  assume  that  they  form 
a  complete  set.  The  two  conditions  can  be  ex¬ 
pressed  in  terms  of  the  eigenfunctions  t Lm  as 

j  =  5LL.6„.m  (15) 

and 

r  tn'i'I*)  =  w-  -  4  W»  -  ID  ■  (16> 

Lm  P  jl'l>V 

To  calculate  two-time  correlations  we  shall 
need  the  Green’s  function  0 (/, ,  I2,  t  j/?,  4>  t0)  which 
is  also  the  conditional  probability  for  the  intensi¬ 
ties  to  be  characterized  by  the  values  4  and  4  at 
time  t  given  their  values  f?,  [°2  at  t0.  To  this  end 
we  note  that  (N (/2 ,  /2,  /  ( /^,  4,  t0)  is  the  solution 
PUl,/.i,t)  of  the  Fokker- Planck  equation  (7)  with 
the  initial  condition  P(4,4,f0)  =  5(/,  -  l?)5(/2  -  4). 

It  follows  immediately  that 

G(4,/„f  ,q,l°2,tB) 

Lm  •  ° 


The  two-time  joint  probaoilitv  density 
/>2(/, , /2, /;/(,/', /')  for  the  intensities  at  two  differ¬ 
ent  tunes  can  be  obtained  easily  now.  Noting  that 
the  stationary- state  probabilities  are  independent 
of  the  origin  of  time  we  obtain 

/>*(/, ,  4,  /  +  T;  = <;(/, ,  4,  r o)ps(i[,i'2) 

=72.'  lmVu  4 If  !,(/,' ,4 ; )c  r , 

L  "i 

F^o.  (18) 

For  discussing  second-order  intensity  correla¬ 
tions,  we  need  the  two-time  joint  probability  densi¬ 
ty,  which  we  have  been  able  to  express  in  terms 
of  certain  eigenvalues  and  eigenfunctions.  For 
discussing  higher-order  intensity  correlations, 
higher-order  probabilities  will  be  required,  which 
also  can  be  expressed  in  terms  of  the  same  set  of 
eigenfunctions  and  eigenvalues.  Our  problem  now 
is  to  solve  for  the  eigenvalues  and  the  eigenfunc¬ 
tions. 

111.  SOLUTIONS  OF  THE  EIGENV  ALUE  PROBLEM 

The  eigenvalue  equation  (12)  can  be  written  in 
a  simpler  form  in  terms  of  two  new  variables  de¬ 
fined  by 

k=4+4,  u  0  |19a) 

JhzJA,  -1*(  -1  (19b) 

Vt*/a) 


We  also  have 
, ,  ,,  u  .  , 


rf/,rf 4  -  2  ■ 

Then  g(m  obeys  the  following  equation: 


((d^.3gcy+„(drn-;<(i  » (i_,-)-i| <1<s0,  (2i) 

Idu  Jtt  v  4,-1-  2.1  ’  )  \  4.4‘:  2  A2  /  4  V  .1  /A  or'  or y Lm 


The  form  of  Eq.  (21)  suggests  a  solution  of  the 
form 

*'«(«.  '•>=«  t*U‘)SLir).  (22 

Substitution  of  this  into  Eq.  (21)  yields  the  allow¬ 
ing  two  uncoupled  equations  for  ltLm{u)  and  Sf(r): 


r 

UC(A-C)  , 

'(.l  -C)‘  3fft  \ 

[  4  a* " 

'  2,1'  "  +' 

4 a*  2  ,1*7 

r 


ft 
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and 


(23b) 


where  ,iL  is  the  separation  constant.  With  the 
choice  ,i  L  -  (L  4)[L  +  2)  =/(/  + 1),  f=0,l,2,  Eq.  (23b) 
is  just  the  Legendre  equation  and  i>\  can  be  ex¬ 
pressed  in  terms  of  P,,  the  Legendre  polynomial 
of  order  l  as 


The  equation  satisfied  by  .  is 


HL  t  2)  t 


-  0  . 


(20 


(241 


The  eigenvalue  equation  (23)  can  be  reduced  to  an 
even  simpler  form  by  making  a  change  of  variable 


u-y *,  - ilu  -  r\lr  , 

K,Jr*)  =  r-^cJr) 


(25) 


This  equation  is  similar  to  the  equation  encountered 
in  the  solution  of  single- mode4  and  two- mode  ring 
laser  problems  ■  except  for  the  /'-dependent  factor 
multiplying  the  eigenvalue.  Unfortunately  Eq. 

(26)  cannot  be  solv'd  analytically  and  numerical 
methods  have  to  be  used.  Finally,  in  terms  of  the 
new  functions  the  two-time  joint  probability  has  the 
form 


x^l+^n'j  S L(i  ')<'‘‘C'"r  ,  T  0. 


(27) 


Higher-order  joint  probabilities  can  be  expressed  similarly.  With  the  knowledge  of  the  two-time  joint 
probability  function,  we  are  now  in  a  position  to  discuss  the  intensity  correlations. 

IV  INTENSITY  CORRELATIONS 

We  shall  consider  only  the  second-order  or  the  two-time  intensity  correlations.  Higher-order  correla¬ 
tions  can  be  discussed  similarly  with  the  knowledge  of  the  higher-order  joint  probabilities.  From  the 
definition  and  by  symmetry,  the  two  autocorrelation  functions  are  equal,  i.e., 


(/,(/)/,  (f  +  D)  =A/e(/)/2(/  +  T)> 

= /// I  nUPiV^v‘  +  T\K,l'z,t)dlxJIjil[dl±. 

The  intensity  cross-correlation  function  is  defined  by 

</2(/)4(t  +  r  r>  =(/,(/'/,(/ +  r» 

=///  J/;/2  />*(/,  ,/„  /  +  r  j/,\  i^Ddi^i'M  ■ 

We  can  evaluate  them  all  together  by  writing 

</,(/>/,</  +  TO  =  ////  +  T\i[,r2,t)di,di2di[di'2 ,  j,y  =  i,2. 

With  the  help  of  Eq.  (27)  and  Eqs.  (18)  and  (19)  we  obtain 

<«,(/)/,.(/  +  D)  =  «  Y.  |jf  <WP00(l+®  <)Rirn\ 

x(~-)1/2e-*i.mr  J'1rf,.|l-(_lV,.jPI(r)  j' d, .'(1  -  (- 1  )'V ] /»,(i-' 
From  the  orthogonality  property  of  the  Legendre  “~ivnrm-2l»,  we  have 
f  '  (1-  (-l)'rjP1(rWr  =  25,0-  ;(-l)'6„  . 


(28a) 


(28b) 


(29) 


),  )./  =  !. 2. 


(30) 


(31) 
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Using  Eqs.  (25)  and  (31)  in  Eq.  (30)  we  obtain  the  following  expression  lor  the  intensity  correlations: 


,/,(/)/,.(/  +  D> 


E 


=  1,2. 


(321 


The  r  integrals  have  to  be  evaluated  numerically. 
The  intensity  correlations,  thus,  can  be  expressed 
as  a  series  of  falling  exponentials  with  coefficients 
which  can  be  evaluated  in  terms  of  the  eigenfunc¬ 
tions  of  Eq.  (27).  The  behavior  of  the  normalized 
intensity  correlation  functions 


*IAT)=ii‘{t)rly,+  T~  ~  1  ■  >’>'  =  1’2  <33) 


is  illustrated  in  Figs.  1  and  2.  Numerical  calcula¬ 
tions  show  that  near  threshold  many  terms  in  the 
series  of  exponentials  contribute  and  a  single  ex¬ 
ponential  approximation  is  rather  poor.  However, 
for  large  excitations  they  are  well  approximated 
by  a  single  exponential.  In  fact. 


r)  <3«, 


is  a  good  approximation  for  .1-'  1.01C.  It  follows 
from  Eq.  (34)  that  with  increasing  excitation  in¬ 
tensity  correlation  time  may  increase  which  is 


f - - 

r Lea riv  reflected  in  Figs.  1  and  2.  These  conclu¬ 
sions  are  similar  to  those  reached  by  Tehram 
and  Mandel  and  Hioe'  for  a  two- mode  laser.  This 
similarity  again  points  out  to  the  fact  that  they 
both  correspond  to  a  system  of  two  neutrally 
coupled  modes  as  noted  in  Ref.  1.  In  view  of  these 
remarks  we  expect  that  the  general  time-dependent 
problem  when  the  losses  of  the  two  modes  are  not 
equal  can  be  discussed  using  a  perturbative  ap¬ 
proach  for  small  differences  in  the  losses  along 
the  lines  of  Ref.  5  and  similar  conclusions  may 
be  reached.  We  shall  not  consider  this  problem 
any  further,  however. 


V  SUMMARY 

We  have  considered  and  solved  the  problem  of 
two-mode  laser  oscillations  on  coupled  transitions 
in  a  A  configuration  of  a  three-level  atomic  sys¬ 
tem.  Two-time  intensity  correlations  have  been 
studied  and  the  expressions  for  them  as  a  series 


FIG.  1.  Normalized  autocorrelation  function  uuif)  as 
a  function  of  T  for  three  operating  points  (A  -  C)/A 
-0.0,  0.007,  0.1.  The  parameters  for  this  figure  are 
A  =  1ms"1  ,  A/B- 10*  and  C  was  varied 


T  'Noturoi  Units) 


FIG.  2.  Normalized  cross-correlation  function  u,j(b 
as  a  function  of  T  for  three  operating  points.  For  para¬ 
meters  see  Fig.  1. 
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of  decreasing  exponentials  have  been  derived.  It 
is  found  that  the  correlation  time  increases  with 
increasing  excitation  and  that  for  even  moderately 
high  excitations  all  second-order  intensity  correla¬ 
tions  are  well  approximated  by  a  single  exponen¬ 
tial. 
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I  he  stejdv  state  i|uantiini  statistical  properties  ot  lieht  in  a  'wo-mode  two-photon  laser  are  presented.  It  in  shown  that 
the  second-order  coherence  functions  violate  a  classical  mequahtv . 


In  many  systems  involving  the  Interaction  between 
light  ami  a  medium,  the  quantum  statistical  proper¬ 
ties  of  light  are  predicted  to  exhibit  some  nonelassieal 
effects  1 1 1.  Hie  most  well-known  example  of  these 
effects  is  photon  antibimcliing  which  was  predicted 
and  then  observed  in  resonance  fluorescence  [2.o| 
and  theoretically  predicted  in  a  multiphoton  absorp¬ 
tion  process  [4.x  | ,  second-harmonic  generation  |o|. 
degenerate  parametric  amplification  [7|.  and  tree- 
electron  laser  [N| .  etc.  We  define  the  degree  of  sec¬ 
ond-order  coherence  of  light  to  be 


<;\r 


(iij  aj  a jOj)- 
(a'j  (ij> 


(1) 


where  a  -  and  a,  arc  the  creation  and  destruction  oper¬ 
ators  ot  the  field  in  the  inode  i  and  the  angular 
brackets  denote  the  ensemble  average.  Hie  photon 
antibunchine  is  exhibited  by  the  fields  which  satisfy 
the  inequality 


(r'.:  1  <  I 


C) 


i.e.  the  degree  of  second-order  coherence  is  less  than 
unity.  I  Ins  is  true  of  the  field  distributions  m  which 
photons  arc  anticorrclatcd. 

In  the  case  of  intensity  measurements  on  the  two 


itcwanh  supported  In  the  l  .S.  Air  I  one  Oltiee  ol 
Seientitie  Keeeareh  (  \l  OSK). 


beams,  we  obtain  | for  a  non  negative  two-mode 
Glauber's  coherent  state  representation  Flu t .  a:  ij 

J  fj  fill]  -  liq  i:!«y|-l: 

X  P(u  | .  Ill  )/,(U| ,  U-1 )  d-i/|  d-l/s  d-t>  ]  d-1'2  >  0  . 

This  leads  to  the  following  quantum  analogue  of  die 
Cauchy's  inequality  |1| 

<;\\'<  .  (?) 

The  violation  of  this  inequality  would  be  expected  in 
systems  where  the  correlation  between  the  two  beams 
is  larger  than  the  correlation  between  the  photons  ot 
the  same  beam.  Tins  nonelassieal  correlation  between 
the  light  beams  lias  been  observed  in  two-photon  cas¬ 
cade  emission  |l)  | . 

In  this  paper  we  give  another  example  where  the 
violation  of  the  inequality  (.'I  could  K’  obseived. 

We  consider  a  two-mode  two-photon  laser  in 
which  lasing  action  is  achieved  by  stimulated  emis¬ 
sion  of  two  photons  in  a  single  atomic  decay  .  The 
laser  model  consists  of  a  coupled  s\  stem  ot  field  and 
identical  two-level  atoms.  Hie  lasinc  levels  ol  the 
atoms  are  assumed  to  have  the  same  pants  under  the 
usual  dipole  approximation  tor  the  two-photon  tran¬ 
sition  to  lake  place. 

Flic  quantum  statistical  properties  ol  the  optical 
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field  ol  a  two-photon  laser  have  been  studied  in  de¬ 
tails  in  recent  years  ( 10  1  7|.  These  include  studies 
on  single-mode  and  two-mode  lasers.  Many  authors 
restrict  themselves  to  the  model  of  a  two-photon  laser 
in  which  the  end  mirrors  are  assumed  to  he  transmit¬ 
ting  at  the  sum  frequency'  ( 10.1 7,14.1  7|.  The  restric¬ 
tion  to  a  two-photon  loss  mechanism  is  made  in  order 
to  retain  the  property  of  detailed  balance.  In  this  pa¬ 
per  we  consider  a  more  realistic  model  of  a  two-mode 
two-photon  laser  in  which  the  cavity  loss  mechanisms 
of  the  two  modes  are  simulated  by  the  single  photon 
processes. 

We  consider  two  field  modes  in  a  resonant  cavity 
interacting  with  a  system  of  identical  two-level  atoms. 
The  atom  field  interaction  may  be  described  by  the 
effective  hamiltonian  in  the  interaction  picture 

//,  =rng(a*'<}|<7i  +  a  d.d,)  .  (4) 

where  g  is  the  coupling  constant  and  the  atom  opera¬ 
tors  a*  and  a"  are  defined  by 

/0  K  /  0  0\ 

a  =  I  .  o~  =  .  (5) 

Vo  o'  Vi  o' 

The  master  equation  for  the  photon  distribution  func¬ 
tion  pin  , .  n  i 1  can  be  obtained  by  a  straightforward 
extension  of  the  Scully- Lamb  theory  [18]  of  a  single 
photon  laser.  The  resulting  equation  is 

d/r(ii,,i<s)  -i4(M|  +  l)(Hi  +  1) 

d7  ~ =  f+  {B/A)ut\V\)0i2  +T) ,,(Wl ’ 

An  ,if  s 

+  ;  .  .  r* — ■  iHn,  -  l.i/s  -  1) 

1  +  (R/A  lit  ,iis  1 

+  (’,(ll,  +  11/1(11,  +  1,11s)-  C|l/i/>(ll, .  Its) 

+  C,(iis  +  1  )/l(U|  .  i/s  +  I )  -  Csiis/'(n, .  ii, 1 ,  ((>) 

where  A  and  R  are  the  gain  and  saturation  parameters 
for  the  two-photon  gam  mechanism  and  C,  (i  =  1 .21 
is  the  loss  parameter  lor  the  fill  mode  which  is  related 
to  the  mode  frequency  co,  and  <J  of  the  cavity  by  the 
relation 

Cj  =  ojj/Q  .  (7) 

This  equation  is  difficult  to  solve  even  in  steady  state 
because  each  element  /i(n, .  n, )  is  coupled  to  /i (n , 

1 .  ii  ,  1 1.  /Mu  ,  +  1  .it, )  and  /)(  ii  | ,  ii,  +  It  which 

represent  two-photon  emission  (one  in  each  mode). 


single-photon  absorption  in  mode  I  and  single-photon 
absorption  in  mode  2  respectively. 

We  can  however  determine  the  second-order  corre¬ 
lation  functions  in  steads  state  when  the  laser  is  oper¬ 
ating  high  above  threshold.  In  tins  legutie  u ,  u  I 
I.  When  this  condition  is  satisfied,  we  can  neglect 
unity  in  comparison  with  the  terms  proportional  to 
R/A  in  the  denominators  of  the  first  two  terms  on  the 
rhs  of  eq.  (ts).  We  then  obtain  (in  steady  state i 

(.4  -/B)/i(ii, .  it,)  +  (.1- '7))/’(ii ,  l.iis  it 
tC,(H,  +  11/1(11,  *•  l.li,)  C,w  ,/’(«, .  II  s  I  IM 

+  f,(n,  + 1  >/i(ii | s  + 1 1  r_s«s/.’(H | .  ii ,  i  -  o . 

Various  moments  of  n ,  and  n  ,.  e.g.  OiVii  s>(g.  d  =  0. 

1 . ...»  can  be  obtained  by  multiplying  eq.  ( M  by 
nail's  and  summing  over  n ,  and  it,.  It  can  be  easily 
shown  that  ( for  f  =  1.2) 

<i!,>  =  A  -;RCi  .  (4) 

</l‘>  =  </!,->“  +  <«,•)  ,  (10) 

(;l]Hs>  =  <«1Xlls)  [1  +  ((ll,>  +  (lls»“  1  ]  .  (11) 

On  recalling  the  definitions  of  Gjy'  [eq.  ( 1  )|.  we  ob¬ 
tain 

C',i)  =  G'sV  =  1  ,  =  1  +  [0i ,>  +  <«,>]  !  .  |12) 

It  is  evident  that  the  Cauchy's  inequality  (.,)  is  violat¬ 
ed  in  the  present  case.  This  nonclassieal  effect  is  how¬ 
ever  very  small,  of  the  order  of(vii,l  +  in,')  '. 

It  is  a  pleasure  to  thank  Professor  M.O.  Scully  for 
his  support. 
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